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PREFACE 

This  technical  report  comprises  publications  made  between  September 
1974  and  the  present.  During  this  period  the  program  was  expanded  to 
include  the  study  of  the  nonlinear  acoustics  of  solids.  Because  of  the 
broadening  of  the  range  of  topics,  in  February  1978  we  decided  to  separate 
out  a single  topic — Ultrasonic  Wave  Reflection  at  Liquid-Solid  Interfaces — 
and  present  it  as  Technical  Report  No.  15.  The  current  technical  report 
overlaps  Technical  Report  No.  15  in  time,  but  not  in  subject  matter.  The 
current  technical  report  gives  our  contributions  to  subjects  other  than 
reflection  of  ultrasonic  waves  at  interfaces.  For  convenience,  the 
report  is  divided  into  three  parts. 

Part  I is  made  up  of  two  papers  dealing  with  critical  mixing  of 
binary  liquid  mixtures.  These  two  papers  essentially  complete  our 
contribution  to  this  subject  for  the  present. 

Part  II  covers  two  different  measurement  techniques.  The  first 
paper  is  a presentation  of  a comparison  of  ultrasonic  pressure  amplitudes 
measured  by  three  different  optical  techniques  and  by  a thermocouple 
probe  with  calculated  amplitudes.  The  measurements  were  made 
simultaneously  in  the  same  ultrasonic  field  in  1958,  but  previously 
were  unpublished.  A revival  of  interest  in  calibration  at  the  Miami 
meeting  of  the  Acoustical  Society  of  America  persuaded  the  author  and 
Floyd  Dunn  that  our  results  would  be  of  interest  to  our  colleagues  and 
thus  prompted  the  publication.  The  second  paper  describes  the  improve- 
ment in  measurement  of  ultrasonic  wave  velocities  in  solids  which  results 
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from  the  use  of  capacitive  transducer.  This  technique  eliminates  the 
perennial  bond  problem. 

Part  III  comprises  six  contributions  to  the  nonlinear  acoustics 
of  solids.  The  first  is  a short  summary  paper  describing  the  temperature 
dependence  of  the  nonlinearity  parameters  of  germanium  and  of  fused 
silica  down  to  77°K.  The  second  is  an  invited  paper,  given  at  the  6th 
International  Symposium  on  Nonlinear  Acoustics  in  Moscow,  which  describes 
the  measurement  technique  and  present  results  on  germanium  down  to  3°K. 

The  third  is  a more  extensive  description  of  both  the  nonlinearity 
parameter  and  the  third-order  elastic  constants  of  germanium  between 
room  temperature  and  3°K.  The  fourth  is  a short  paper  given  at  the  ICA 
in  Madrid  on  the  nonlinearity  parameter  of  four  different  types  of  fused 
silica.  The  fifth  is  a more  extensive  account  of  the  results  of  our 
experiments  on  the  nonlinearity  of  fused  silica.  The  sixth  and  final 
paper  is  a summary  paper,  given  at  the  International  Symposium  on  Nonlinear 
Acoustics  in  Paris,  which  gives  the  temperature  dependence  of  the 
nonlinearity  parameters  and  third-order  elastic  constants  of  germanium 
and  copper  and  makes  an  attempt  to  relate  them  to  other  quantities  of 
significance  in  solid  state  physics. 
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Ultrasonic  study  of  critical  mixing  of  hexane  and 
0,  |3'-dichloroethyl  ether* 

Elizabeth  M.  Bains  and  M.  A.  Breazeale 

Department  of  Physics,  The  University  of  Tennessee,  Knoxville.  Tennessee  37916 
(Received  8 February  1974) 

The  velocity  and  attenuation  of  ultrasonic  waves  in  the  critical  temperature  region  of  a binary 
mixture  of  hexane  and  J3,  jS’-dichloroethyl  ether  have  bcjn  measured  at  1,  3,  i,  7,  and  9 MHz.  By 
assuming  that  the  attenuation  per  wavelength  is  made  up  of  two  terms,  aX  =Af  ,I,  + Bf.  we  obtain 
good  agreement  with  the  Fixman  theory  for  both  aX  and  a//2.  The  magnitude  of  B,  however,  is 
only  partially  accounted  for  by  our  measured  values  of  shear  viscosity. 
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INTRODUCTION 


Examination  of  the  attenuation  maximum  for  ultrasonic 
waves  in  the  critical  mixing  region  of  binary  liquids  was 
begun  by  Lucas  and  Biquard,  1 who  attempted  to  explain 
the  observed  maximum  by  analogy  with  the  scattering  of 
visible  light  from  density  fluctuations.  This  explanation 
had  been  successful  for  critical  opalescence.  Scattering 
of  ultrasonic  waves  in  the  critical  mixing  region  was  in- 
vestigated by  Brown  and  Richardson. 2 Later  Brown2  re- 
ported that  scattering  accounts  for  the  major  part  of  the 
exc  jss  attenuation  away  from  the  critical  point  but  a 
much  smaller  part  as  the  critical  point  is  approached. 

A relaxation  due  to  rearrangement  of  molecular 
clusters  was  proposed  by  Chynoweth  and  Schneider,  4 but 
Anantaraman,  Walters,  Edmonds,  and  Pings2  and 
D’Arrlgo,  Mistura,  and  Tartaglia8  reported  that  a single 
relaxation  does  not  account  for  the  frequency  dependence 
of  the  data. 


The  attenuation  of  sound  in  fluids  is  classically  as- 
sumed to  depend  upon  two  transport  properties,  shear 
viscosity  and  thermal  conductivity: 


aX  = 


ITU)  /4 
c*Po\3 


V + 


(y-  1 )K\ 

C , ) ’ 


(1) 


where  u>  is  the  angular  frequency,  c is  the  sound  veloc- 
ity, Po  is  the  undisturbed  density,  t)  Is  the  viscosity  co- 
efficient, y is  the  ratio  of  specific  heats,  K is  the  ther- 
mal conductivity  coefficient,  and  Cp  is  the  specific  heat 
at  constant  pressure.  If  the  thermal  conductivity  term 
is  negligible,  then  the  attenuation  per  wavelength  is  di- 
rectly proportional  to  the  viscosity  coefficient,  and  is 
also  a linear  function  of  frequency.  Although  an  in- 
crease in  viscosity  at  the  critical  mixing  temperature 
was  measured  early  in  the  study  of  binary  liquids, 7 and 
a theoretical  explanation  of  it  has  been  offered, 8 the  ef- 
fect of  this  Increase  heretofore  has  been  assumed  to  be 
negligible  in  explanations  of  the  increase  in  attenuation 
in  the  critical  region.  The  magnitude  of  this  increase 
in  viscosity  is  typically  5%-40%,  whereas  the  magnitude 
of  the  Increase  in  ax  is  typically  factors  of  2- 10.  In  the 
data  to  be  presented,  however,  the  effect  of  shear  vis- 
cosity is  noticeable. 

In  addition  to  this  shear  viscosity,  there  is  evidently 
an  Increase  in  bulk  viscosity  in  the  critical  region.  The 
effect  of  the  bulk  viscosity  on  the  attenuation  of  sound 


has  been  considered  by  Kawasaki  and  Tanaka, 9 who 
show  that  bulk  viscosity  gives  an  attenuation  per  wave- 
length' which  has  the  same  form  as  that  derived  by  Fix- 
man.  10 

The  theory  derived  by  Fixman10  using  an  excess  heat 
capacity  rather  than  bulk  viscosity  as  the  mechanism 
will  be  used  in  the  interpretation  of  our  results.  This 
theory  gives  an  /"1/4  frequency  dependence  for  aX,  the 
temperature  dependence  of  aX,  and  the  temperature  and 
frequency  dependence  of  the  velocity  in  the  critical  re- 
gion. 

Fixman  derives  the  excess  attenuation  and  dispersion 
of  sound  from  the  interaction  between  concentration  fluc- 
tuations near  the  critical  point  and  a sound  wave  passing 
through  the  region.  The  temperature  variations  asso- 
ciated with  the  adiabatic  sound  wave  alter  the  concentra- 
tion fluctuations,  and  these  fluctuations  decay  toward  the 
equilibrium  distribution  according  to  the  diffusion  equa- 
tion. This  decay  distorts  the  temperature  variation  of 
the  sound  wave,  and  causes  absorption  of  energy  from 
the  wave. 

In  recent  years  there  have  been  a number  of  experi- 
mental tests  of  the  Fixman  theory.  In  a number  of 
cases8’8'11"12  it  was  observed  that  the  frequency  depen- 
dence of  a//2  of  the  data  followed  that  predicted  by  Fix- 
man.  10  Anantaraman  et  al.  * reported  that  although  a//* 
showed  the  frequency  dependence  predicted  by  Fixman, 
ax  did  not. 

The  purpose  of  the  present  work  is  to  show  how  the 
frequency  dependence  of  both  aX  and  a//1  of  critical 
mixture  of  hexane  and  0,  /3'-dlchloroethyl  ether  (com- 
monly called  chlorex)  can  be  brought  into  agreement 
with  the  Fixman  theory.  This  is  done  by  assuming  an 
excess  attenuation,  whose  frequency  dependence  is  the 
same  as  that  caused  by  viscosity,  in  addition  to  the  ex- 
cess attenuation  resulting  from  the  Fixman  mechanism. 

THEORY 

Fixman  calculates  an  excess  heat  capacity  per  unit 
volume  due  to  variations  in  entropy,  expressed  in  terms  of 
Fourier  components  of  the  radial  distribution  function.  The 
perturbation  of  the  equilibrium  components  caused  by  tem- 
perature variations  in  the  sound  wave  is  evaluated  by  in- 
troducing into  the  diffusion  equation  a heat  capacity  that 
is  dependent  both  on  local  temperature  and  on  Fourier 
components  of  the  radial  distribution  function  that  are 
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exponential  functions  of  time.  The  excess  heat  capacity 
is  obtained  in  terms  of  a complex  integral  and  is 


where 

F(tl)=  d~  i,2j“  u>*(wi  + 1 H-uT’  + ieV2*  l)Prf«\  (3) 
d-  kz(Ii/w)'/2  , (4) 


dependence  of  the  experimental  aX.  According  to  equa- 
tions (10)  and  (11),  the  excess  attenuation  per  wavelength 
ax  should  depend  on  f'vi.  As  will  be  shown,  the  experi- 
mental data  do  show  such  a dependence;  however,  it  is 
necessary  first  to  remove  the  linear  frequency  depen- 
dence to  be  expected  from  shear  viscosity.  In  order  to 
do  this,  it  is  necessary  to  introduce  a plot  of  a//2  as  a 
function  of  frequency.  This  allows  one  to  separate  out 
the  shear  viscosity  contribution  as  a constant  term. 


and 

kjJ^  »|lj(N|*«1) 
w,  2m  pc  2 ' 

The  integration  variable  in  Eq.  (2)  is  w = k/n,  and  the 
parameters  are  Boltzmann’s  constant  kB,  temperature 
T,  critical  temperature  Tc,  parameters  of  the  radial 
distribution  function  k and  a,  angular  frequency  w,  dif- 
fusion constant*,  and  molecular  mass,  molecular  vol- 
ume, volume  fraction,  number  density,  and  mass  frac- 
tion, /«,,  4 >,,  «2  and  c , , Cj,  respectively,  for 

liquids  1 and  2. 

The  imaginary  part  of  the  integral  in  Eq.  (3)  deter- 
mines the  sound  absorption  and  the  real  part  of  that  in- 
tegral determines  the  dispersion  of  the  speed  of  sound. 

The  calculation  of  the  attenuation  and  dispersion  of 
sound  is  made  by  assuming  that  CP  and  Cv  are  increased 
by  the  same  amount  by  the  composition  fluctuations 


CP  = C°P  + 


CY 


= C»  + 

«l  + «2 


(6) 


where  N0  is  Avogadro’s  number  and  Cp  and  C°  are  heat 
capacities  in  the  absence  of  critical  composition  fluctu- 
ations. 


The  sound  wave  propagating  along  the  x axis  is  given 
by 

ip  - exp i(u>x/cc  - u >/)  , (7) 

where  <■„  is  the  complex  speed  of  sound.  It  is  related  to 
c0,  the  speed  of  sound  In  the  absence  of  critical  com- 
position fluctuations,  by 


cc*c0 


(8) 


The  real  part  of  cc  is  the  speed  of  sound  c.  The  imag- 
inary part  of  iA>/ce  is  the  attenuation  o.  The  attenuation 
per  wavelength  is  2vcaa/w.  Using  Eq.  (2),  the  speed  of 
sound  and  the  attenuation  per  wavelength  can  be  written 
in  terms  of  the  complex  integral  Fid)  as 


c = c*{l-  i WRe[/'(d)]} 


(9) 


and 


(10) 


where 


4 71*  (*i  + nt)  \UJ; 


(id 


and  R is  the  gas  constant. 

The  most  direct  comparison  of  the  Flxman  theory  with 
experiment  comes  from  a determination  of  the  frequency 


EXPERIMENTAL  TECHNIQUES 

A critical  mixture  of  chlorex  and  hexane  was  sealed 
in  a cylindrical  glass  tank.  A drawing  of  the  apparatus 
is  shown  in  Fig.  1.  Measurements  of  the  ultrasonic  at- 
tenuation and  velocity  were  made  by  single  crystal  pulse 
techniques,  with  the  pulses  propagating  between  a piezo- 
electric quartz  crystal  in  the  transducer  and  the  air- 
backed  mica  sheet  in  the  reflector.  The  transducer 
could  be  translated  through  the  liquids  relative  to  the  re- 
flector, and  the  distance  through  which  it  moved  could 
be  measured  with  a 3-in.  dial  gauge  to  an  accuracy  of 
±0.001  in. 

The  electronic  circuits  used  to  generate  and  analyze 
the  ultrasonic  pulse  are  shown  in  Fig.  2.  An  rf  pulse 
generated  by  an  Arenberg  Model  PG650-G  pulsed  os- 
cillator was  fed  into  both  an  echo  detecting  circuit  and  a 
frequency  measuring  circuit.  In  the  echo  detecting  cir- 
cuit a diode  switching  circuit  isolated  the  pulsed  oscilla- 
tor from  the  crystal  and  also  protected  the  detector  by 
clipping  the  initial  pulse  while  minimizing  attenuation  of 
the  low  amplitude  echoes.  The  1 MHz  signals  were  fed 
directly  to  an  oscilloscope  (Tektronix  Type  555),  but  the 
higher  frequency  signals  were  detected  by  a Hewlett- 
Packard  855  Spectrum  Analyzer  before  going  to  the  os- 
cilloscope. For  the  attenuation  measurements  the  height 
of  the  first  echo  was  measured  with  a voltage  comparator 
unit  (Type  W)  for  different  transducer  positions. 
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FIG.  2.  Electronic  system. 


For  the  velocity  measurements,  echoes  from  double 
pulses  generated  by  the  Arenberg  oscillator  were  over- 
lapped to  measure  wavelength  by  a modified  interferom- 
eter technique.  The  frequency  counting  circuit  used  a 
mixer  (HP10504A/B)  to  mix  the  short  (i.e. , 20  psec)  rf 
pulses  with  a cw  signal  generated  by  an  HP8601A  Signal 
Generator.  The  difference  frequency  output  signal  from 
the  mixer  was  displayed  with  the  oscilloscope  and  the 


frequency  of  the  generator  was  adjusted  for  zero  beat. 
The  frequency  of  the  generator  was  then  counted  by  a 
General  Radio  1191-13  Counter. 

Temperature  regulation  was  provided  for  the  liquids 
under  investigation  through  heat  exchange  liquids  sur- 
rounding the  tank.  The  temperature  of  the  heat  exchange 
liquids  was  regulated  by  aHaake  Circulator.  This  regula- 
tion was  ±0.  02  °C  inside  the  circulator.  The  temperature 
stability  of  the  binary  liquids  was  much  better  than  this 
because  they  were  inside  a carefully  insulated  container. 
The  temperature  of  the  binary  liquids  was  measured  by 
the  glass-encapsulated  iron-constantan  thermocouple 
shown  in  Fig.  1 and  displayed  on  a strip  chart  recorder. 

After  assembly,  the  apparatus  was  tested  by  making 
a measurement  of  a//2  for  distilled  water  between  0 and 
23.  5 °C.  The  measured  values  of  a//2  fell  within  the 
experimental  scatter  of  the  data  of  Smith  and  Beyer.14 
They  stated  an  average  probable  error  of  5%  for  their 
measurements. 

In  the  calibration  data,  as  well  as  in  the  critical  mix- 
ing data,  a diffraction  correction  was  made. 15  This  cor- 
rection was  as  large  as  50%  at  1 MHz,  and  was  only  1% 
at  9 MHz.  The  intermediate  frequencies  had  a correc- 
tion between  these  extremes. 

Velocity  of  ultrasonic  waves  in  distilled  water  mea- 
sured between  0 and  20  °C  exhibited  a mean  deviation  of 
0. 17%  from  the  curve  given  by  Wilson. 16 

Viscosity  measurements  were  made  using  an  Ostwald 
capillary  viscometer  in  a bath  whose  temperature  was 
controlled  by  an  Artronix  Model  5301  Temperature  Con- 
troller. The  viscometer  was  calibrated  by  using  dis- 
tilled water  as  a reference  liquid. 
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FIG.  3.  Ultrasonic  wave  velocity  In  the  single-phase  region  above  the  critical  temperature  12.7*C. 


J.  Cham.  Phyi.,  Vol.  61,  No.  3,  1 August  1974 


L.  IVI.  001113  ailU  IVI.  /-\.  uiuaicaic,  uiuojuim,  itioui  niiAtuiuu 


I /.H  I 


FIG.  4.  Ultrasonic  wave  velocity  in  chlorex  as  a function  of 
temperature. 

RESULTS  AND  DISCUSSION 
Velocity 

Measured  values  of  the  ultrasonic  wave  velocity  in  the 
critical  temperature  region  are  shown  in  Figs.  3(a)-3(e) 
for  the  frequencies  1,3,5,  7,  and  9 MHz,  It  is  seen  that 
the  velocity  in  the  single-phase  temperature  region  in- 
creases as  the  critical  temperature  is  approached.  This 
behavior  is  the  reverse  of  that  predicted  by  the  Flxman 
theory,  but  is  reconcilable  with  it, 8 Both  chlorex  and 
hexane  have  a negative  temperature  coefficient  larger 
than  that  for  the  mixture  at  the  critical  temperature. 

The  temperature  coefficient  of  hexane  is17  - 4.4  m/ 
sec  °C.  The  velocity  in  chlorex  as  a function  of  temper- 
ature is  given  in  Fig.  4,  where  the  temperature  coeffi- 
cient is  toand  to  be  - 3.  57  m/aec  °C.  (The  data  points 
for  5 MHz  are  shown  in  the  figure;  those  for  1 and  9 
MHz  are  not.  The  scatter  at  these  frequencies  was  sim- 
ilar to  that  at  5 MHz. ) The  temperature  coefficient  of 


the  velocities  in  Fig.  3 is  ~ - 4 m/sec  "C  away  from  the 
critical  point.  It  decreases  to  ~ - 5 m/sec  “C  at  the 
critical  point.  Thus,  it  is  seen  that  the  temperature  co- 
efficient of  the  critical  mixture  becomes  smaller  than 
that  of  the  noncritical  mixture  (or  of  the  separated  liq- 
uids) in  agreement  with  the  Flxman  theory. 

Attenuation 

In  Figs.  5(a)-5(e)  are  presented  measured  values  of 
the  attenuation  per  wavelength  aA  as  a function  of  tem- 
perature in  the  single-phase  temperature  region.  The 
solid  curves  are  a best  fit  of  the  data  points,  assuming 
the  functional  form  of  the  Fixman  theory  given  in  Eq. 

(10)  with  an  additional  term  to  account  for  the  noncritical 
attenuation: 

aA  = tfWIm[F(rf)]+  const.  (12) 

The  tables  of  Kendig  et  al. 18  were  used  to  evaluate  the 
integral  Fid).  The  values  obtained  by  extrapolating  these 
curves  to  the  critical  temperature  12.  7 °C  are  shown  in 
Table  I as  aA^,,,^,.  This  procedure  depends  upon  a 
tacit  assumption  that  all  mechanisms  contributing  to  the 
critical  attenuation  (shear  viscosity,  etc. ) have  approxi- 
mately the  same  temperature  dependence.  The  fit  of  the 
experimental  data  is  taken  as  an  indication  of  the  validity 
of  this  approximation;  however,  an  improvement  in  the 
approximation  is  possible. 

Measurement  of  the  shear  viscosity  as  a function  of 
temperature  allows  one  to  Introduce  the  temperature 
dependence  of  the  noncritical  attenuation  through  use  of 
Eq.  (1).  (The  temperature  dependence  of  velocity  is  ob- 
tained from  Fig.  3. ) Both  viscosity  and  velocity  data 
were  obtained  from  linear  extrapolation  of  the  asymp- 
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FIG.  5.  Attenuation  per  wavelength  In  the  single-phase  region  above  the  critical  temperature  12.7"C.  Solid  line  Is  Flxman  theory. 
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TABLE  I.  Parameters  used  In 


According  to  Eq.  (10),  a A.. Is  proportional  to  f'Vi. 

The  remaining  attenuation  per  wavelength  has  been  la- 
beled “vtscltke”  because  It  Is  proportional  to  /,  as  Is  the 
viscosity  contribution. 

These  values  of  ax  can  be  expressed  In  terms  of  a/f* 
by  dividing  by  c/.  Values  of  (a//1),,,*,,  are  given  In 
Table  I.  Since  viscosity  contributions  to  a/f*  are  con- 
stant In  frequency,  one  can  plot  a/f*  as  a function  of 
frequency  and  Identify  the  vtscosttyllke  contribution  from 
an  evaluation  of  the  Intercept.  We  chose  to  use  the  fre- 
quency dependence  of  f'vt  expected  from  the  Flxman 
contribution  as  abscissa  In  Fig.  7.  The  fact  that  the 
data  fall  on  a straight  line  with  a nonzero  Intercept  In- 
dicates that  a/f*m,m  Is  adequately  expressed  by  the  sum 
of  two  terms: 

W/*)»w  = A’r*/i+B’.  (15) 

The  coefficient  A'  Is  Identified  with  the  Flxman  contri- 
bution; B’  Is  viscositylike.  Since  the  magnitude  of  Bf 
Is  Independent  of  frequency,  this  Intercept  value  can  be 
subtracted  at  each  frequency  to  Isolate  the  contribution 
of  the  Flxman  mechanism. 

The  effect  on  aAT,.T„  of  subtracting  the  vlscosltyllke 
contribution  Is  Indicated  In  Fig.  8 In  which  we  have 
plotted  aX,ntm  both  before  and  after  subtracting  the  vls- 
costtylike  contribution  obtained  from  Fig.  7.  As  can  be 
seen,  after  correction  the  data  follow  a better  straight- 
line  behavior  as  a function  of  f‘u*,  In  agreement  with 
the  Flxman  theory.  The  uncorrected  data  show  a cur- 
vature not  seen  In  aXriiaM.  Use  of  any  value  of  B'  other 
than  the  Intercept  value  of  Fig.  7 Introduces  a similar 
curvature  In  aA.. 

One  can  attempt  to  explain  this  vlscosltyllke  contri- 
bution by  assuming  that  all  viscous  mechanisms  In- 
crease In  proportion  to  shear  viscosity  In  the  critical 
region.  The  viscosity  Increased  by  6%  In  the  critical 
region,  as  shown  In  Fig.  8.  This  would  lead  to  an  ex- 
cess attenuation  Indicated  by  the  cross-hatched  region  In 


totlc  values  away  from  the  critical  point.  For  example, 
the  assumed  noncrltical  contribution  to  viscosity  Is  In- 
dicated by  a dotted  line  In  Fig.  6,  In  which  we  have 
plotted  experimental  data  for  the  chlorex- hexane  mix- 
ture. 

It  one  considers  only  shear  viscosity,  the  attenuation 
per  wavelength  calculated  from  Eq.  (1)  is  much  smaller 
than  that  observed.  Since  we  are  unable  to  calculate  the 
bulk  viscosity  contribution,  we  simply  multiply  by  the 
factor  required  to  bring  the  viscous  attenuation  to  80% 
of  the  attenuation  calculated  from  Eq.  (12)  at  30  °C. 

(The  trend  of  the  data  Indicates  thai  the  noncrltical  con- 
tribution to  aX  has  approximately  this  value  at  that  tem- 
perature. The  use  of  100%  Instead  of  80%  alters  the 
magnitudes  of  calculated  coefficients  by  a small  amount, 
but  does  not  alter  the  conclusion. ) Using  the  extrapolated 
noncrltical  temperature  dependence  of  v and  c,  the  non- 
crltical attenuation  Is  extrapolated  to  12. 7 °C.  The  cal- 
culated curves  are  superimposed  on  Fig.  5.  The  critical 
contribution  from  shear  viscosity  and  velocity  would 
raise  these  curves  by  approximately  5%  at  the  critical 
temperature.  On  the  scale  of  Fig.  5,  this  difference  Is 
lass  than  the  width  of  ths  Uns  labelled  a\  noncrltical. 

The  difference  between  the  two  curves  in  Fig.  5 at 
12. 7 *C  has  been  defined  as  the  excess  attenuation  per 
wavelength  resulting  from  critical  phenomena  as  follows: 

almwa  " aWttl»l  • (*3) 


TNs  excess  attenuation,  given  In  column  3 of  Table  I 
Is  observed  to  be  made  up  of  two  parts  which  have  dif- 
ferent frequency  dependence: 

* a*rtlBM+a*VlM,n*  • 


f**  (MHrr** 

FIG.  7.  Critical  point  values  of  the  exoess  or//1  fitted  to  f*'* 
dependence. 


FIG.  S.  Sheer  viscosity  of  the  binary  mixture  shove  the  orltloal 
temperature  12.7‘C. 
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FIG.  8.  Critical  point  values 
of  the  attenuation  per  wave- 
length fitted  to  the  frequency 
dependence  predicted  by  Fix- 
man. 


wavelength  to  bulk  viscosity  effects;  however,  such  an 
assumption  would  be  inconsistent  with  the  results  of 
Kawasaki  and  Tanaka, a who  show  that  In  the  critical  re- 
gion the  bulk  viscosity  produces  an  aX  which  has  the 
same  form  as  the  Ftxman  integral  [Eq.  (3)). 

Note  added  in  Proof:  A recent  publication  by  Peter 
H.  Rogers  and  A.  L.  VanBuren  [j.  Acoust.  Soc.  Am. 
68,  724  (1074)|  gives  a simple  closed-form  expression 
which  makes  the  correction  of  attenuation  data  for  the 
effects  of  diffraction  quite  direct.  In  addition,  this  ex- 
pression is  valid  over  a wider  range  of  distances  than 
the  graphical  correction  of  Seki,  Granato,  and  Truell 
(Ref.  IS). 


Fig.  7,  which  Is  less  than  20%  of  the  viscositylike  con- 
tribution. We  see,  therefore,  that  simple  assumptions 
about  the  viscous  mechanisms  do  not  allow  us  to  account 
for  the  excess  attenuation.  But  whatever  its  origin,  the 
remaining  excess  attenuation  appears  to  have  a frequency 
dependence  which  Is  the  same  as  that  of  viscosity. 

CONCLUSION 

Measured  values  of  velocity  and  attenuation  of  ultra- 
sonic waves  In  the  critical  temperature  region  of  a 
chlorex- hexane  mixture  have  been  compared  with  the 
theory  of  Flxman.  The  velocity  data  show  a trend  which 
Is  consistent  with  the  Flxman  predictions.  The  attenua- 
tion data  are  found  to  be  described  by  a sum  of  two 
terms.  The  frequency  dependence  of  one  term  Is  In 
agreement  with  that  predicted  by  the  Flxman  theory. 

The  frequency  dependence  of  the  second  term  Is  that  ex- 
pected from  shear  viscosity;  however,  the  magnitude  of 
the  viscositylike  term  Is  greater  than  can  reasonably  be 
expected  from  the  critical  Increase  In  shear  viscosity 
for  this  mixture.  Although  there  should  be  a similar 
frequency  dependence  resulting  from  thermal  conductiv- 
ity, this  contribution  Is  expected  to  become  vanishingly 
small  In  the  critical  region. 9 It  would  be  tempting  at 
this  point  to  attribute  the  remaining  attenuation  per 
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Recently  we  published  measurements  of  the  attenua- 
tion and  velocity  of  ultrasonic  waves  In  the  critical 
temperature  region  of  the  binary  system  hexane  and 
P,  ^'-dlchloroethyl  ether. 1 By  assuming  a "viscosity- 
like”  contribution  to  the  critical  attenuation  In  addition 
to  the  noncrltlcal  attenuation  (which  also  has  the  fre- 
quency dependence  expected  for  attenuation  due  to  vis- 
cosity), we  were  able  to  obtain  agreement  with  the  the- 
ory of  Fixman  for  the  critical  aA  as  well  as  a//*.  Pre- 
vious authors,  in  fitting  data  to  the  Fixman  theory, 
found  agreement  with  a//*,  but  not  with  aX,*~*  In 
this  note  we  demonstrate  that  the  assumption  of  a 
critical  region  "vlscosity-llke"  contribution  in  the  data 
for  several  other  mixtures  identifies  a term  in  the  crit- 
ical aX  which  agrees  with  the  / ',/4  frequency  dependence 
of  the  Fixman1  theory. 

Figure  1 gives  a//1  as  a function  of  f "*'*  for  the 
binary  mixtures  listed.  These  data  were  taken  from 
Refs.  1-6.  In  every  case  the  agreement  between  data 
points  and  the  straight  line  shown  (the  least  squares  fit 
of  the  data)  is  satisfactory,  confirming  the  agreement  of 
a/f*  with  the  Fixman  theory  reported  for  these  data. 


The  nonzero  Intercept  value  Isolates  a part  of  a/f*  that 
is  Independent  of  frequency,  as  is  the  attenuation  due  to 
For  all  the  data  except  that  of  Ref.  I,  the 
intercepts  Include  the  noncrlttcal  attenuation  as  well  as 
the  critical  region  "viscosity-like”  contributions.  The 
noncrltlcal  attenuation  was  subtracted  from  each  data 
point  in  Ref.  1 so  that  the  Intercept  is  only  the  “viscos- 
ity-like” contribution. 

In  Fig.  2 are  given  the  same  data  plotted  in  the  form 
<*A  as  a function  of  /'1/4,  which  is  the  frequency  depen- 
dence predicted  by  Fixman.  The  values  of  aA  were 
calculated  by  subtracting  the  Intercept  values  shown  in 
F i.».  1 from  the  appropriate  experimental  values  of 
a/f*  and  multiplying  by  cf.  As  can  be  seen,  the  data 
points  in  Fig.  2 agree  quite  well  with  the  straight  lines, 
with  the  possible  exception  of  the  data  of  D’Arrlgo, 
Mlstura,  and  Tartaglla.  * Several  of  their  data  points 
are  also  seen  not  to  lie  exactly  on  the  line  of  Fig.  1,  A 
smaller  Intercept  on  the  a/f  * curve  would  lead  to  a 
more  linear  aA  curve;  however,  for  consistency  we 
chose  to  define  the  Intercept  by  the  least  squares  fit  of 
a/f * as  a function  of  f*1*,  The  curvature  of  the  aA 
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FIG.  1.  Measured  frequency  dependence  of  a//2  for  binary 
liquid  mixtures,  (a)  Reference  1 . (Chlorex  is  the  common 
name  for  (3,  0'  -dlehloroethyl  ether.)  (b)  Reference  2.  (c)  Ref- 
erence 3.  Id)  Reference  4.  (e)  Reference  5.  (f)  Reference  (i. 
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curves  was  found  to  be  very  sensitive  to  the  magnitude 
of  intercept  used. 

In  Table  I we  compare  the  value  of  intercept  found 
from  Fig.  1 (for  each  set  of  data  other  than  that  of 
Ref.  1)  with  the  value  for  noncritical  attenuation  ob- 
tained for  those  data  by  the  experimenter.  For  the 
data  of  Ref.  1,  as  stated  above,  the  figure  shows  val- 
ues of  excess  a//2  with  the  noncritical  contribution 
subtracted,  so  the  existence  of  an  intercept  in  that  case 
clearly  indicates  a critical  region  “viscosity-like”  term. 
Other  experimenters  expressed  their  total  a//2  in  the 
form 

*/ = Af"*'*  ImfF(d)]  + B(T),  (1) 

where  Im|F(r/)l  is  the  imaginary  part  of  the  integral 
given  by  Fixman,7  d depends  on  the  temperature,  and 
A and  B(T ) are  constants  obtained  by  fit  of  the  data  to 
the  temperature  dependence  predicted  by  Fixman.  The 
term  B{T)  is  identified  as  the  noncritical  part  of  a//2. 
D'Arrigo,  Mlstura,  and  Tartaglia  gave  the  value  of 
B(T)  listed  in  Table  I from  their  fit;  the  other  values 
are  average  values  for  all  frequencies  calculated  from 
experimental  values  of  a//2  and  the  empirical  A.  In 


FIG.  2.  Frequency  dependence  of  oA  for  binary  liquid  mix- 
tures. (a)  Reference  1.  (b)  Reference  2.  (c)  Reference  3.  (d) 

Reference  4.  (e)  Reference  5.  (Velocity  was  assumed  to  be 

1250  m/sec.  Another  value  would  change  the  slope  of  the  line 
slightly,  but  would  make  no  change  In  the  scatter  of  the  data. ) 

(f)  Reference  (i. 


each  case  the  intercept  (which  is  the  value  required  to 
isolate  the  Fixman  term  in  aA)  is  larger  than  the  pre- 
viously calculated  noncritical  value  of  a//2.  This 
would  have  remained  true  even  if  we  had  used  a smaller 
intercept  to  bring  the  data  of  D’Arrigo,  Ml6tura,  and 
Tartaglia  into  bette*-  alignment. 

In  summary,  the  procedure  we  have  outlined  isolates 
that  part  of  the  critical  attenuation  that  is  accounted  for 
by  the  theory  of  Fixman.  For  the  cases  given  above,  it 
also  shows  that  a “viscosity-like”  term  exists  in  the 
critical  region.  For  the  hexane  and  0,  ^-dichloroethyl 
ether  system,  it  was  found1  that  the  critical  region  in- 
crease in  shear  viscosity  was  not  sufficiently  large  to 
account  for  the  magnitude  of  the  “viscosity-like”  term. 
Whether  or  not  this  is  true  for  the  other  cases  cannot  be 
answered  at  present  because  of  the  lack  of  data  on  ex- 
cess critical  point  viscosity  for  the  other  systems. 


TABLE  I.  Comparison  between  Intercept  values  of  or//2  from 
f Ig.  1 and  the  total  noncritical  o//2  obtained  by  the  experi- 
menter. 


Binary  mixture 

Intercept 

(aecVcm) 

fl(T)(=  a//2  noncritical) 
(iec2/cm) 

Nitrobenzene -n-hexane* 

148x  10*IT 

126x10’" 

Nltrobenzene-iao-octaneb 

435x10-" 

315x10’" 

Tr>thyl  amine-water* 

565X10"" 

343x10’" 

n-  Heptane— o-tol  uidl  ned 

348  x 10"" 

186X10’" 

Aniline -cyclohexane* 

148x10’" 

100X  10’" 

•Reference  2. 

•Reference  3. 

•Reference  4. 

•Reference  5. 
•Reference  6. 

•Research  supported  In  part  by  the  U.  S.  Office  of  Naval  Re- 
search. 
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Comparison  of  methods  for  absolute  calibration  of  ultrasonic 
fields* 

M.  A.  Breazeale 

Physics  Department,  University  of  Tennessee,  Knoxville.  Tennessee  3 79] 6 

F.  Dunn 

Bioacoustics  Research  Laboratory.  University  of  Illinois.  Urbana,  Illinois  61801 
(Received  10  December  1973) 

Methods  of  determining  the  absolute  ultrasonic  pressure  amplitude  under  identical  circumstances  are 
compared  Methods  used  are  (1)  radiation  force  on  a small  sphere,  (2)  thermoelectric  probe,  and  (3) 
three  optical  techniques. 

Subject  Classification:  35.80,  80.30 


At  the  recent  Society  meeting  in  Miami  (28  November- 
1 December  1972),  the  Technical  Committee  on  Physical 
Acoustics  sponsored  three  sessions  on  ultrasonic  energy 
and  power  measurements.  These  sessions  included 
state-of-the-art  invited  papers,  tutorial  papers,  round 
table  discussions,  and  contributed  papers  on  current  in- 
vestigations. The  attendance  and  discussion  participa- 
tion clearly  showed  that  such  measurements,  and  their 
interpretation,  are  of  fundamental  importance  in  a va- 
riety of  ultrasonic  research  programs  During  the  round 
table  discussions  and  discussions  following  formal  pres- 
entations, a question  was  raised  about  the  degree  to 
which  various  absolute  measurement  procedures  com- 
pare. Though  considerable  interest  was  exhibited  in 
this  point,  no  published  works  were  recalled,  extempo- 
raneously, for  citation.  The  present  authors  were  re- 
minded that  they  carried  out  such  a comparison  involv- 
ing five  different  methods,  at  two  laboratory  sites,  in 
1958,  but  that  the  results  had  not  been  published.  We 
report  herein  our  findings  in  the  hope  that  we  may  an- 
swer some  of  the  questions  raised  in  Miami. 


FIG.  1.  Calibration  of  a 990-kHz  quartz  transducer:  (A)  cal- 
culated output  with  «n-4.77xl04i  (B)  thermocouple  probei  (C) 
optical  linage  broadening  (refraction);  CD)  light  diffraction;  CE) 
decrease  of  iigbl  Intensity  (refraction). 
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The  circumstances  of  this  comparison  involved  the 
research  programs  at  Michigan  State  University  under 
the  direction  of  the  late  E.  A.  Hiedemann  and  at  the 
University  of  Illinois  under  the  direction  of  the  late  W. 

J.  Fry.  The  former  group  had  been  concerned  with 
three  optical  methods,  viz.,  the  refraction  (image broad- 
ening) method, 1 the  diffraction  method, * and  the  refrac- 
tion (decrease  of  light  intensity)  method, 3 while  the  lat- 
ter group  had  been  employing  the  method  of  radiation 
force  on  a small  solid  sphere4  and  the  transient  thermo- 
electric method.5  While  each  group  had  found  favorable 
comparison  among  the  methods  they  embraced,  more 
extensive  treatment  was  considered  desirable.  The  pro- 
cedure involved  the  Illinois  group  calibrating  one  of  its 
thermocouple  probes  against  the  radiation  force  detector 
(small  solid  sphere)  and  transporting  it  to  the  Michigan 
State  campus  where  the  former’s  experimental  arrange- 
ment was  duplicated.  Thus  the  thermocouple  probe  was 
operated  by  the  Michigan  State  group  in  the  same  way  as 
in  Illinois  and  its  output  compared  with  that  detected  by 
the  three  optical  methods  mentioned  above.  Figure  1 
illustrates  this  comparison  in  water  for  the  output  from 
a 6-cm-diam  quartz  transducer,  fundamental  thickness 
resonant  frequency  990  kHz.  Curve  A is  the  output  cal- 
culated under  the  assumption  that  the  transducer  is  a 
piston- like  vibrator.  The  piezoelectric  constant  en 
= 4.77x10*  was  used  in  this  calculation.  Curve  B is  the 
average  of  values  from  the  thermocouple  probe,  alter 
correction  for  temperature  variation  c i the  output.  Curve 
C is  from  the  optical  image  broadening  method.1  Curve 
D is  from  the  light  diffraction  method.8  Curve  E is  from 
the  decrease  of  light  intensity  method.*  It  is  seen  that 
the  experimental  results  exhibit  a total  range  of  approx- 
imately ± 10%  about  the  mean,  and  that  this  mean  is  ap- 
proximately 27%  below  that  calculated  from  the  voltage 
applied  to  the  transducer.  Because  of  the  assumptions 
made  in  the  calculation,  one  should  consider  this  cal- 
culated value  to  be  the  maximum  value  possible  under 
ideal  conditions.  The  true  value  should  be  less  than  this, 
as  the  experimental  results  show. 

Although  further  refinements  of  the  optical  methods* 
have  reduced  the  scatter  among  the  data,  and  the  avail- 
ability of  lasers  now  makes  the  optical  methods  more 
direct,  we  present  these  results  as  a guide  in  estlmat- 
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Elimination  of  transducer  bond  corrections  in  accurate 
ultrasonic-wave  velocity  measurements  by  use  of  capacitive 
transducers* 

John  H.  Cantrell,  Jr  and  M.  A.  Breazeale 

Department  o / Phyiuv  I' he  Untventty  of  Tenntnee.  Knoxville.  Tennessee  37916 
i Received  4 May  1^76.  revised  5 October  1976) 

A ..apautive  driver -capacitive  detector  system  for  generation  and  detection  of  ultrasonic  waves  has  been 
developed  I hi*  eliminate*  the  necessity  of  bonding  piezoelectric  transducers  to  solid  samples  With  the 
capacitive-dnver  capacitive-detector  system,  free-free  boundary  conditions  exist  at  the  sample  surfaces 
and  longitudinal  ultrasonic-wave  velocities  in  solids  can  be  measured  accurately  without  correcting  for 
ultrasonic  wave  phase  shift*  due  to  sample-bonded  transducer  interfaces  The  capacitive  driver  has  a mica 

1 dielectric  which  increases  the  breakdown  potential,  hut  maintains  the  free-free  boundary  conditions  at  the 

solid  specimen  surfaces  This  allows  for  a largcr-amplitude  ultrasonic  signal  to  be  generated  in  the  sample 
than  is  possible  with  an  air-gap  capacitive  driver.  This  improves  the  precision  of  the  measurement.  The 
accuracy  of  the  method  is  comparable  with  that  of  bonded-transducer  methods,  after  bond  corrections  are 
made. 

PACS  numbers  43.33. Yb.  43.33.Cg 


INTRODUCTION 

Ultrasonic  pulse-echo  techniques  are  widely  used  in 
the  measurement  >A  elastic -wave  velocities.  Three 
fundamental  pulse -echo  methods  for  measuring  phase 
velocities  in  solids  have  been  reported.  They  are  the 
gated  double-pulse  superposition  method  of  Williams 
and  L«imb.‘  the  pulse-superposition  method  of  Mc- 
Skimin,1  and  the  echo-overlap  technique  of  May1  and 
Papadakls. 4 

Until  recently,  it  has  been  necessary  to  generate  the 
ultrasonic  pulses  and  to  detect  the  resulting  echoes  by 
applying  bursts  of  rf  to  a piezoelectric  transducer 
bonded  to  the  specimen  surface.  Such  an  arrangement 
necessitates  making  a correction  in  the  velocity  mea- 
surement due  to  ultrasonic-wave  phase  shifts  at  the 
specimen-transducer  interface.  In  some  situations, 
the  effect  of  this  coupling  between  specimen  and  trans- 
ducer leads  to  considerable  error  If  neglected.1 

In  order  to  eliminate  the  bond  problem  in  compres- 
sional-wave  measurements,  the  air-gap  capacitive  de- 
tector previously  used  to  measure  amplitudes  of  ultra- 
sonic waves*'7  was  modified  to  function  as  an  air-gap 
capacitive  driver  as  well. ' The  air-gap  capacltlve- 
driver-capacitive-detector  combination  allowed  the 
sample  ends  to  vibrate  with  free-free  boundary  con- 
ditions exactly.  No  bond  corrections  were  necessary 
because  there  were  no  bonds. 

It  is  the  purpose  of  this  paper  to  report  on  the  develop- 
ment of  a dielectric  capacitive  driver  to  replace  the 
air-gap  capacitive  driver.  This  new  driver  is  capable 
of  generating  longitudinal  ultrasonic  pulses  whose  am- 
plitudes are  much  larger  than  those  in  the  previous 
system.  We  show  that  free-free  boundary  conditions 
still  exist  at  the  sample  surface  with  this  driver  and 
that  the  error  in  the  velocity  measurements,  without 
making  any  corrections.  Is  comparable  to  the  error  In 
velocity  measurements  made  with  bonded  transducers, 
after  corrections  for  bond  phase  shifts. 


I.  THE  CAPACITIVE  DRIVER-CAPACITIVE  DETECTOR 
SYSTEM 

Figure  1 shows  the  mechanical  parts  of  the  dielec- 
tric-capacitive -drive  r-air-gap-capacitive -detector 
system  with  the  sample,  in  place.  The  sample  rests  on 
the  grounded  outer  portion  of  the  detector  assembly. 

The  grounded  outer  portion  of  the  dielectric  capacitive 
detector  is  recessed  approximately  10  pm  so  that  the 
electrode  and  sample  face  form  a parallel -plate  capaci- 
tor. The  detector  electrode  is  kept  at  a dc  bias  of  the 
order  of  150  V.  The  ultrasonic  vibration  of  the  sample 
face  varies  the  gap  spacing  (and  hence  the  capacitance), 
thereby  generating  an  alternating  electrical  signal  be- 
tween the  electrode  and  ground.  This  signal  is  am- 
plified and  displayed  on  the  oscilloscope.  The  capaci- 
tive detector  has  been  refined  to  the  point  that  it  i9 
capable  of  giving  absolute  displacement  amplitudes  for 
30-MHz  ultrasonic  waves  having  amplitudes  as  small 
as  l<r4  A.’ 


FIG.  1.  Mechanical  parts  of  the  capsclllve-driver- 
capacttlve-detector  system. 
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FIG.  2.  Dielectric  driver 
arrangement. 


( t ) SAMPLE 


In  the  air-gup  capacitive  driver  the  inner  electrode 
is  also  recessed  approximately  10  pm  to  form  a paral- 
lel plate  capacitor  with  the  sample  surlace,  but  in  the 
dielectric  capacitive  driver  a 7.  5-pm-thick  piece  of 
mica  <3)  is  sandwiched  between  the  sample  surface 
[outer  electrode  < 1 ) j and  the  inner  electrode  (6)  as 
shown  in  Fig.  2.  Between  the  electrodes  and  the  cen- 
tral mica  dielectric  (3)  two  air  layers  ({2)  and  (41) 
are  present  because  of  slight  irregularities  of  the  con- 
tact surfaces. 

The  capacitive  driver  is  not  dc  biased,  but  a pulsed 
sinusoidally  varying  rf  voltage,  V * F0slnu>/,  applied 
between  the  electrodes  j (1 ) and  (5)J,  produces  sinu- 
soidally varying  electric  fields  Eu  F.z,  and  £,.  It  has 
been  demonstrated10  that  the  generation  of  ultrasonic 
waves  in  this  situation  is  due  to  the  effect  of  electro- 
static forces  acting  directly  on  the  ultrasonic-wave 
propagation  medium.  Hence,  if  E,  is  the  electric  field 
near  the  surface  of  the  sample,  the  electrical  force  per 
unit  area  />  at  this  surface  is 

. (1) 

where  «0  is  the  dielectric  permittivity  of  air.  Since 
£|  * £,0sinwf , (2) 

then 


/>*  }c0£,l0stnIu>/  =»  Uo^lo  0 - cos2w/). 


(3) 


For  present  purposes,  It  is  adequate  to  ignore  the  con- 
stant term  In  Eq.  (3)  and  consider  only  the  cos2u>/  term. 
One  sees  that  the  pressure  on  the  sample  surface  is 
applied  at  (nice  the  frequency  of  the  rf  voltage.  Thus, 
the  ultrasonic  wave  generated  in  the  sample  has  twice 
the  frequency  of  the  applied  rf  voltage. 

II.  the  experimental  system 

A.  Method  of  measurement 

The  velocity  measurements  were  made  with  the  equip- 
ment arrangement  as  shown  in  Fig.  3.  The  gated 
double-pulse  superposition  method  of  Williams  and 
Lamb1  was  employed.  A pulsed  ultrasonic  signal  is 
fed  Into  the  sample  from  the  capacitive  driver  followed 
bv  another  pulsed  ultrasonic  signal  phase  locked  but 
delayed  In  lime  with  respect  to  the  first.  The  delay  is 
adjusted  such  that  superposition  of  the  desired  echoes 
from  the  two  pulse  trains  is  achieved.  The  resulting 
signal  is  received  at  the  other  end  of  the  sample  by 
the  capacitive  detector,  passes  through  a wide-hand 
i.f.  amplifier  whose  center  frequency  is  20  MHz,  and 
is  displayed  on  an  oscilloscope.  Utilizing  the  fact  that 
a continuous  change  in  ultrasonic  frequency  results  In 


a succession  of  superimposed  echo  maxima  and  minima, 
the  ultrasonic  velocity  is  calculated  from 

v * 2wf(A//Atf),  (4) 

where  is  the  change  in  ultrasonic  frequency  corre- 
sponding to  An  minima,  / is  the  sample  length,  and 
w is  an  integer  (2  in  our  experiments)  giving  the  rela- 
tive position  of  the  initial  pulse  train  and  the  delayed 
pulse  train.  To  enhance  the  sensitivity  In  detecting 
the  superimposed  echo  minima,  a boxcar  integrator 
was  incorporated  into  the  system  between  the  i.  f. 
amplifier  oulpul  and  the  oscilloscope. 

B.  A tuned  transmission-line  voltage  transformer 

The  capacitive  transducer  is  driven  with  a pulsed  rf 
signal.  The  frequency  is  varied  around  15  MHz  to  pro- 
duce successive  resonances.  To  obtain  the  voltages 
necessary  for  adequate  ultrasonic  amplitudes,  a select 
length  of  transmission  line  was  used  to  act  as  a tuned 
voltage  transformer  between  the  gated  amplifier  and 
the  capacitive  driver. 

If  .S'  is  the  length  of  the  transmission  line,  f Is  the  rf 
dnve  frequency,  Vs  is  the  voltage  across  a purely  re- 
active load  of  magnitude  \ZL\  (capacitive  driver),  V0 
is  the  voltage  out  of  the  gated  amplifier,  and  Z0  is  the 
characteristic  impedance  of  the  transmission  line, 
then  fora  lossless  transmission  line 


V , » Fosin0  csc(«f>  - wS/r), 

where 

0 - -tan"1  \Zt\/Z0  , 
u»«2 »/  , 


(5) 

(6) 
(7) 

and  c is  the  velocity  of  propagation  of  electromagnetic 
waves  in  the  transmission  line.  Thus,  for  the  proper 
lengths',  Fs/F0  becomes  very  large.  But  the  impedance 
looking  out  of  the  gated  amplifier  into  the  transmission 
line  behaves  as 

Zu*-jZ0t3in(<t>  -u>S/c).  (8) 

Thus,  for  the  condition  that  V9/V0  becomes  large,  Zu, 
and  consequently,  V9,  become  small.  For  this  reason, 
the  length  of  transmission  line  is  chosen  such  that  Z, 


FIG.  3.  Hlock  diagram  of  experimental  apparatus. 
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TABLE  1.  Sample  lengths  and  comparative  values  of  A//Aw  anti  velocity  for  air-gap  driver  measurements  and 
dielectric  driver  measurements. 


Air  gap  driver 

Dielectric  driver 

Sample  length 

A fl  An 

Velocity 

A// Aw 

Velocity 

Sample 

(cm! 

ulu  - MHz) 

<x  104  cm/sec! 

(*10’2  MHz) 

(*  10-  cm/sec) 

S-I 

1.251)3  *0.0002 

5.8926  to. 0017 

5.9364  *0  0020 

5.8925*0.0006 

5. 9363  *0.0012 

S-II 

1.2586 *0.0002 

5. 0957  * 0.0013 

5.9362  * 0.0016 

5.8958  »0.0005 

5.9363  *0.0011 

SW-I 

1.2500  t 0.0003 

5.0009  *0.0017 

5.9525  t 0.0023 

5.9100  *0.0005 

5.9526  * 0.0016 

SW-II 

1.2587  1 0.0002 

5.9121  : 0.0013 

5.9532  *0.0016 

5.9122  *0.0004 

5.9533  s 0.0011 

GE151 

1.2721  *0.0001 

5.6315*0.0016 

5.9378  t 0.0017 

5.8345  *0.0006 

5.9378  *0.0008 

looks  inductive  for  our  operation  frequency.  A variable 
capacitor  is  then  placed  in  parallel  with  Zu  to  form  a 
parallel  resonant  circuit.  This  arrangement  increases 
the  intensity  of  the  generated  ultrasonic  signal  by  ap- 
proximately 30  dl3  over  that  of  the  untuned  transmission 
line. 

The  primary  limitation  on  the  generated  ultrasonic- 
wave  amplitude  is  set  by  dielectric  breakdown  in  the 
capacitive  transducer.  However,  this  limitation  is  not 
as  serious  as  one  might  assume.  As  the  spacing  of 
the  electrodes  in  air  is  decreased,  the  breakdown  po- 
tential goes  up  almost  exponentially.*  In  the  air-gap 
driver  we  were  able  to  use  fields  as  large  as  200  kV/cm 
without  breakdown.  With  the  mica  dielectric  driver, 
we  have  been  able  to  use  fields  larger  than  1000  kV/cm. 

III.  EXPERIMENT  AND  RESULTS 

In  order  to  determine  the  effect  of  contact  between 
the  mica  dielectric  and  the  sample  surface,  velocity 
measurements  were  made  on  five  fused  silica  samples 
at  room  temperature.  Each  sample  was  measured  ten 
times  with  the  dielectric  capacitive  driver  and  ten 
times  with  the  air-gap  capacitive  driver.  The  results 
of  the  measurements  are  shown  in  Table  I.  The  aver- 
age values  and  standard  deviations  of  A// An  and  lon- 
gitudinal velocities  |Eq.  (4)|,  as  well  as  sample  length 
1 are  given. 

Previous  work  comparing  velocity  measurements  with 
bonded  transducers  and  air-gap  transducers  has  con- 
firmed experimentally  that  the  air-gap  capacitive-d  river- 
capacitive-detector  system  eliminates  the  effects  of 
transducer  bonding.*  The  present  agreement  between 
the  air-gap-  and  dielectric-driver  measurements  of 
Af/ Aw  for  each  sample  confirms  within  experimental 
error  that  the  dielectric  driver  also  eliminates  the 
effects  of  transducer  bonding. 

A.  Analysis  of  error 

The  standard  deviation  in  \f/ A//  with  the  dielectric 
driver  (approximately  0.008'7/)  is  significantly  lower 
than  with  the  air-gap  driver  (approximately  0.025%) 
because  of  the  much  larger  ultrasonic  signal  generated 
by  the  dielectric  driver.  The  null  frequencies  of  the 
superimposed  pulse  trains  were  reproducible  to  within 
t 30  Hz,  which  is  approximately  two  parts  in  10*.  Hie 
relatively  large  standard  deviations  in  the  velocity  cal- 
culations In  Table  I are  due  to  lack  of  parallelism  of 
the  sample  surfaces,  which  is  reflected  in  the  relatively 
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large  uncertainties  in  the  sample  lengths  listed. 

With  electrically  nonconducting  samples,  it  is  neces- 
sary to  coat  the  end  surfaces  with  a thin  film  of  con- 
ducting material  to  function  as  the  ground  plates  of  the 
capacitive  transducers.  We  found  that  approximately 
1000  A of  copper  is  ideal  for  our  purposes.  According 
to  our  calculations,  this  coating  introduces  an  error  of 
no  more  than  0.003%  in  the  velocity  measurements. 

Corrections  in  the  velocity  measurements  due  to  dif- 
fraction were  calculated11”13  and  serve  to  decrease  the 
values  listed  in  Table  I by  approximately  0. 005%. 

It  is  found  that  the  total  uncertainty  in  compressional- 
wave-velocity  measurements  with  the  dielectric  capaci- 
tive-driver- capacitive-detector  system  is  comparable 
to  the  uncertainty  in  velocity  measurements  with  sys- 
tems using  bonded  transducers.  l,*“*,u 

B.  Advantages  of  the  new  system 

The  advantages  of  the  previous  capacitive-drlver- 
capacitlve-detector  system  have  been  enumerated.* 

They  are : 

(1)  The  tedious  calculation  of  the  effect  of  ultrasunic- 
wave  phase  shift  at  the  bonded  transducer-sample  in- 
terface has  been  eliminated,  with  no  sacrifice  in  ex- 
perimental accuracy. 

(2)  The  effect  of  the  transducer  resonance  has  been 
eliminated.  Capacitive  transducers  have  a flat  fre- 
quency response. 

(3)  The  effect  of  a variation  in  electronic  loading  and/ 
or  the  change  in  thickness  of  the  bonded  transducer  due 
to  temperature  changes  In  cryogenic  velocity  measure- 
ments have  been  eliminated.  (These  variations  often 
force  a change  in  rf  drive  frequency  of  the  bonded  trans- 
ducer to  maintain  a given  pulse  shape. ) 

(4)  Accuracy  of  attenuation  measurements  is  im- 
proved. It  is  well  known,  especially  in  measurements 
at  cryogenic  temperatures,  that  the  effect  of  transducer 
bonds  often  sets  a limit  on  the  accuracy  of  attenuation 
measurements  by  introducing  temperature -dependent 
phase  shifts  in  the  received  signals.  Such  effects  are 
eliminated  with  the  capacitive  transducer. 

In  addition  to  these  advantages,  the  introduction  of  a 
dielectric  into  the  capacitive  driver  allows  one  to  work 
with  large-amplitude  signals.  As  we  have  shown,  this 
dielectric  does  not  affect  the  free-free  boundary  con- 
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ditlons  if  one  is  careful  to  keep  it  dry.  Therefore,  the 
accuracy  of  the  previous  capacitive  driver  Is  main- 
tained, while  the  larger-umplitude  signals  allow  one  to 
improve  the  precision  of  the  measurement. 
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ULTRASONIC  NONLINEARITY  PARAMETERS 


Breazeale  M A 


University  of  Tennessee,  Knoxville,  Tennessee, 
USA 


A consistent  definition  of  the  nonlinearity  parameters  of  gases, 


liquids,  and  solids  can  be  made, 
nonlinear  wave  equation 


By  using  an  appropriate  form  of  the 


3 u 

°0  2 

3x 


32u 

K2  77  + 

3x 


the  nonlinearity  parameter  < **  - ( 


(3K2  + K3) 

3K2  + K3 


3u  3 u 
3*  3x2 


) can  be  recognized  as  the 


negative  of  the  ratio  of  the  coefficient  of  the  nonlinear  term  to  that 
of  the  linear  term.  The  nonlinearity  parameter  has  considerable  sig- 
nificance for  solids*»2  as  well  as  for  liquids  and  gases.  For  iso- 
tropic solids  and  the  [100]  direction  in  single  crystals,  K 


Thus , x 


i + W 


-o  ♦ HiU). 

C11 


Cu  and 


For  liquids  < ■ 2 + B/A.  For 


3 vlll* 
gases 

For  example,  we  calculate  room  temperature  values:  for  air,  x 
2.4,;  for  water  x - 7.0;  for  copper  * 5.25;  for  germanium 


-9.2.  In  our  harmonic  dis- 


+6 


. He  1110] 
» Ge  1111] 

. Ge  [100] 


100  200  300 

Temperature  (*K) 


x . _ . * 2.78.  But  for  fused  silica 

tdrtiOn  experiments,  the  negative 
nonlinearity  parameter  implies  that 
the  waveform  approaches  a "backward 
sawtooth"  as  the  wave  progresses. 

In  addition,  the  magnitudes  of  x for 
fused  silica-^  seem  to  be  sensitive 
to  small  amounts  of  impurities.  An 
OH  content  of  1200  ppm  led  to  a 
measured  x - -9.64.  The  same  fused 
silica  with  only  5 ppm  OH  had  a 
measured  x - -9.25. 

Figure  1 shows  recently  measured 
nonlinearity  parameters  of  ger- 
manium^ and  fused  silica  plotted  as 
a function  of  temperature.  Consis- 
tent with  other  single  crystals 
measured  to  date^,  germanium  x's  are 
almost  independent  of  temperature. 

That  for  the  fused  silica  sample  is 
not.  [Research  supported  in  part  by 
the  U.S.  Office  of  Naval  Research.] 
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Acoustical  Nonlincarities  in  Solids 
As  a Function  of  Temperature 

by 

M.  A.  Breazcale 
Department  of  Physics 
The  University  of  Tennessee 
Knoxville,  Tennessee  37919 


ABSTRACT 


Measurement  of  the  fundamental  and  of  the  second  harmonic  of 
an  initially  sinusoidal  ultrasonic  wave  in  a solid  allows  one  to  compute 
I its  nonlinearity  parameter.  From  this  nonlinearity  parameter  one  can 

isolate  simple  combinations  of  third  order  clastic  constants. 

We  have  developed  apparatus  for  making  measurement  of  the 
fundamental  and  of  the  second  harmonic  of  an  initially  sinusoidal  30  MHz 
ultrasonic  wave  between  300° K and  3°K. 

This  apparatus  and  current  results  for  copper,  germanium  and  fused 
silica  are  discussed.  It  is  shown  that  the  finite  amplitude  distortion 

It  S 

in  fused  silica  is  such  that  at  the  discontinuity  distance  one  would 
have  repeated  rarefaction  shocks.  Both  copper  and  germanium  produce 
distortion  which  would  lead  to  repeated  comprcssional  shocks. 

In  contrast  to  an  assumption  usually  made  in  testing  theories  of 
the  solid  state,  the  third  order  elastic  constants  can  change  as  a 
function  of  tcmpcraturc--as  much  as  150»  between  300°K  and  3°K  in 
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A detailed  investigation  of  the  subject  of  nonlinear  acoustics  is  very 
important  to  the  study  of  all  acoustical  phenomena.  This  is  because  technological 
progress  has  made  possible  the  generation  of  large  amplitude  acoustical  signals, 
and  these  signals  of  necessity  propagate  through  a nonlinear  medium.  There  are 
only  a very  few  media  which  are  truly  linear,  and  even  those  which  appear  to  be 
linear  exhibit  this  behavior  only  under  very  special  circumstances.  An  exact 
statement  of  what  is  meant  by  the  words  "linear  medium"  is  surprisingly  elusive, 
but  this  problem  need  not  concern  us  at  the  present  time.  Rather,  I want  to 
discuss  for  a moment  some  very  profound  implications  of  the  study  of  nonlinear 
acoustics  to  the  subject  of  solid  state  physics. 

The  derivation  of  the  nonlinear  wave  equation  to  describe  acoustical 
nonlinearities  in  solids  can  be  quite  complicated- for  the  general  case.  However, 
for  special  crystal  symmetries  the  equations  take  a much  simpler  form.  The 
nonlinear  equations  for  Isotropic  solids  and  cubic  crystals  are  as  simple  as 
that  for  fluids.  As  a matter  of  fact,  these  equations  are  identical  to  the 
equation  for  fluids  within  the  approximation  involved  in  making  a Taylor 
expansion  of  the  nonlinear  equation  for  fluids  and  dropping  higher-order  terms. 
This  is  the  reason  one  can  speak  meaningfully  about  nonlinearity  parameters  of 
solids  as  well  as  nonlinearity  parameters  of  fluids. 

It  is  through  these  nonlinearity  parameters  that  one  can  make  a fundamental 
connection  between  nonlinear  acoustics  and  solid  state  physics.  As  a matter  of 
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fact,  nonlinear  acoustics  offers  the  possibility  to  measure  quantities  in 
solid  state  physics  which  cannot  be  measured  by  other  techniques.  Of  course, 
in  order  for  the  acoustical  quantities  to  be  useful  to  people  in  solid  state 
physics,  a consistent  definition  of  them  must  be  made,  and  this  has  been  done 
by  Thurston.^- 

In  my  laboratory  we  have  devoted  considerable  effort  to  the  development  of 
some  new  techniques  for  measuring  the  nonlinearity  parameters  of  solids. 

: 

Although  initially  it  was  sufficient  to  be  able  to  measure  nonlinearity 
parameters  at  room  temperature,  it  soon  became  apparent  that  our  results  would 

I 

be  much  more  useful  to  solid  state  physics  if  we  could  measure  as  a function  of 
temperature  down  to  temperatures  near  the  absolute  zero.  The  most  obvious 
reason  for  this  is  that  theories  of  the  solid  state  usually  do  not  take  into 
account  the  effect  of  thermal  motion  on  elastic  properties,  and  thus  are  made 
under  the  tacit  assumption  of  absolute  zero  of  temperature. 

When  one  measures  as  a function  of  temperature,  the  apparatus 
of  necessity  becomes  somewhat  complicated.  But  the  principles  can  remain  simple. 
I should  like  to  describe  the  apparatus  we  have  developed  in  the  past  several 
years — emphasizing  the  principles  involved — and  to  give  some  of  our  recent 
results  on  germanium  as  an  indication  of  the  solid  state  physics  data  available 
from  studies  of  nonlinear  acoustics. 

An  initially  sinusoidal  ultrasonic  wave  becomes  distorted  as  it  propagates 
through  a nonlinear  medium.  Fourier  analysis  of  this  distorted  wave  shows  that 
the  second  harmonic  is  the  largest  overtone.  It  grows  linearly  with  propagation 
distance,  and  is  proportional  to  the  square  of  the  frequency .and  the  square  of 
the  amplitude  of  the  fundamental  component.  In  order  to  write  down  the 


I 


differential  equation  whose  solution  behaves  in  this  fashion,  one  can 
specialize  the  equations  to  crystals  having  cubic  symmetry  and 
write 


P/sU  - K.(U  + 3U  U ) + K,U  U 
0 tt  2 aa  a aa  3 a aa 


where  the  subscripts  denote  partial  differentiation  with  respect  to  the  time 
and  the  space  variable.  The  solution 

U - U(0)  + U(1\  (2! 

under  the  assumption  that  a sinusoidal  wave  is  generated  at  a = 0,  has  the 
proper  behavior.  Measurement  of  the  amplitude  of  the  fundamental  component 

U(0)  - A sin(ka-a)t)  O) 


and  of  the  second  harmonic  component 


3Ko  + K- 

-( — on )(AK)  a cos  2(ka-ut) 


gives  all  the  information  necessary  to  calculate  the  nonlinearity  parameter 


3K2  + K3 


for  a solid  if  one  confines  his  attention  to  the  three  principal  directions  in 
the  cubic  lattice,  or  any  direction  in  an  Isotropic  medium. 

Before  I go  on  to  define  what  K2  and  are,  let  me  digress  for  a moment  and 
point  out  that  the  amplitude  of  the  second  harmonic  (Eq.  4)  has  a minus  sign  in 
front  of  it.  It  turns  out  that  for  most  media  is  negative  and  has  a magnitude 
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greater  than  3^.  This  means  that  the  second  harmonic  is  positive.  We  reported 

2 

in  London  last  year  that  an  exception  to  this  general  rule  exists.  In  fused 

silica  at  room  temperature  was  assumed  to  be  positive,  with  the  result  that 

in  fused  silica  the  minus  sign  is  not  concelled.  This  means  that  the  second 

harmonic  is  generated  180°  out  of  phase  with  that  generated  in  other  solids. 

Since  that  time  we  have  developed  a phase  sensitive  detector  to  measure  the 

relative  phase  of  the  second  harmonic  in  copper  and  in  fused  silica.  The  result 

is  an  unequivocal  proof  that  the  second  harmonics  are  indeed  180°  out  of  phase 
3 

in  the  two.  This  phase  sensitive  detector  is  now  used  routinely  to  monitor 
the  phase  of  the  second  harmonic  as  we  measure  as  a function  of  temperature. 
Therefore,  in  our  data  there  is  no  possible  ambiguity  in  sign  of  the  second 
harmonic. 

Let  us  now  return  to  the  definition  of  and  K^.  For  cubic  crystals 
the  expressions  for  and  are  given  in  Table  I.  As  can  be  seen,  for  the 
principal  directions  and  are  linear  combinations  of  second-order  and  third- 
order  elastic  constants,  respectively.  The  second-order  elastic  constants  can 
be  obtained  from  measurement  of  ultrasonic  wave  velocity.  Thus,  measurement 
of  the  second  harmonic  of  an  initially  sinusoidal  ultrasonic  wave  gives  values 
of  these  combinations  of  third-order  elastic  constants.  For  purposes  of  data 
analysis,  I want  to  point  out  that  the  third-order  elastic  constant  appears 

In  all  of  the  expressions  for  K^.  Thus,  these  are  not  the  simplest  combinations 
of  third-order  elastic  constants  available  from  our  measurements,  since  can 

be  subtracted  from  the  expression  for  in  the  [110]  direction.  A similar 
simplification  can  be  made  In  the  expression  for  In  the  [111]  direction,  as  I 
will  Ghow  shortly. 
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Measurement  of  the  amplitude  of  the  fundamental  and  of  the  second 

harmonic  component  is  made  by  use  of  a special  form  of  capacitive  microphone. 

We  have  measured  the  fundamental,  the  second  harmonic  and  the  third  harmonic 

of  a pulsed  30  MHz  ultrasonic  wave.  A drawing  of  the  capacitive  microphone 

designed  for  measurement  as  a function  of  temperature  is  shown  in  Figure  1. 

An  x-cut  quartz  transducer  generates  the  ultrasonic  wave  which  travels  upward 

in  the  sample.  The  ends  of  the  sample  are  optically  flat  and  parallel.  The  end 

of  the  sample  and  the  detector  electrode  form  a parallel  plate  capacitor,  whose 

spacing  can  be  adjusted  pneumatically  by  changing  the  pressure  differential 

P2~P^.  This  adjustment  is  necessary  because  differential  thermal  expansion 

changes  the  spacing  as  the  temperature  is  changed.  Typically,  this  spacing  is 

adjusted  from  3 to  10  microns  as  the  pressure  differential  is  changed  by  one 

atmosphere.  With  such  small  spacing  of  the  electrodes,  and  an  applied  d.c. 

bias  of  the  order  of  150  volts,  the  microphone  is  very  sensitive.  We  can 

measure  microvolts  of  signal  generated  as  the  pulsed  ultrasonic  wave  impinges 

on  the  end  of  the  sample.  This  means  that  ultrasonic  wave  amplitudes  as  small 
-4 

as  10  Angstroms  can  be  measured,  as  can  be  seen  in  Figure  2,  which  is  a plot 

of  measured  values  of  the  third  harmonic  as  a function  of  the  cube  of  the  funda- 

-3  -4 

mental  amplitude.  The  data  points  are  found  at  0.1  x 10  Angstroms  - 10 
Angstroms.  The  slope  of  this  curve  gave  a value  of  third-order  elastic  constants 
which  agreed  with  values  measured  by  other  techniques  to  within  approximately 
10%.  There  were  no  adjustable  parameters  in  the  measurement. 

A schematic  diagram  of  the  apparatus  for  measuring  as  a function  of 
temperature  is  shown  in  Figure  3.  The  position  of  the  capacitive  microphone 
la  labelled  "detector."  A pulsed  30  MHz  signal  is  generated  at  the  gated 
amplifier.  This  signal  drives  an  x-cut  quartz  transducer  to  generate  an 
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ultrasonic  wave  in  the  sample.  The  capacitive  microphone  on  the  opposite  end 

of  the  sample  produces  an  electrical  signal  proportional  to  the  ultrasonic  ' 

wave  amplitude.  This  signal  is  amplified  and  is  measured  on  an  oscilloscope. 

In  order  to  measure  the  phase  of  the  second  harmonic  we  replaced  the  60 
MHz  amplifier  by  a slightly  more  complicated  frequency  doubler  phase  sensitive 
detector  shown  in  Figure  4.  In  the  upper  area  the  30  MHz  signal  is  fed 

through  a phase  shifter  and  a frequency  doubler.  The  60  MHz  signal  is  then  ' 

caused  to  interact  coherently  in  the  phase  sensitive  detector  with  the 

acoustically  generated  second  harmonic.  The  output  is  then  observed  on  an  ; 

oscilloscope.  As  1 have  indicated,  the  results  were  unequivocal  for  fused  silica 

and  we  are  monitoring  the  phase  as  we  measure  as  a function  of  temperature  in  J 

f 

order  to  detect  any  changes  in  the  sign  of  the  nonlinearity  parameter. 

We  have  measured  the  values  of  for  copper  and  germanium  single  crystals 
down  to  77°K,  as  shown  in  Figure  5.  These  combinations  of  third-order  elastic 
constants  are  linear  functions  of  temperature  over  this  temperature  range.  As 
1 can  be  seen,  some  of  the  lines  have  an  appreciable  slope,  so  that  even  on  the 

basis  of  these  measurements  we  can  state  that  it  is  inadequate  to  assume  that 
the  third-order  elastic  constants  are  independent  of  temperature.  This 
assumption  is  often  made  in  theories  of  the  solid  state. 

Recently  we  have  been  able  to  extend  the  measurements  in  germanium  down 
to  much  lower  temperatures — down  to  2.9°K.  The  results  are  shown  in  Figure  6. 

On  our  first  run  we  were  uncertain  about  the  phase  of  the  second  harmonic,  and 
hence  the  sign  of  in  the  temperature  interval  between  15*K  and  approximately 
50*K.  We  had  observed  the  fact  that  the  Gruneisen  parameter,  which  is  a 
function  of  the  third-order  elastic  constants,  becomes  negative  in  this  region, 
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then  is  positive  again  at  lower  temperatures.  A second  run  with  the  phase 
sensitive  detector  produced  the  solid  data  points  in  the  [110]  direction.  The 
uncertainty  has  been  removed.  These  data  represent  our  measurement  of  these 
combinations  of  third-order  elastic  constants. 

These  combinations  of  third-order  elastic  constants  are  not  the  simplest 
combinations  available  from  our  data,  however.  As  I have  already  indicated, 
and  certain  other  combinations  of  third-order  elastic  constants' appear  in  more 
than  one  of  these  curves,  and  can  therefore  be  subtracted  out.  The  result  is 
shown  in  Figure  7.  These  combinations  of  third-order  elastic  constants  are  the 
simplest  combinations  available  from  our  data  without  making  any  assumptions 
about  the  behavior  of  certain  of  them.  We  see  that  is  a linear  function 

of  temperature  as  far  as  we  can  measure,  but  the  other  combinations  are  strong 
functions  of  temperature  below  approximately  25 °K.  Of  necessity,  the  behavior 
shown  for  the  combination  + ^C144  + ®C456  can  only  if  one  or 

more  of  these  constants  becomes  positive  at  very  low  temperatures.  The  impli- 
cations of  this  statement  to  solid  state  physics  are  not  truly  appreciated  at 
the  present  time.  In  addition,  the  behavior  of  the  Gruneisen  parameter  of 
germanium  is  not  mirrored  in  these  data.  Measurements  of  other  combinations  of 
third-order  elastic  constants  will  have  to  be  made  before  one  can  be  certain 
about  any  single  constant  other  than  The  most  direct  method  for  making 

these  measurements  is  not  clear  at  the  moment.  Nevertheless,  I think  it  is 
correct  to  say  that  Investigations  in  nonlinear  acoustics  are  making  possible 
a much  more  detailed  understanding  of  solid  state  physics. 
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Third-order  elastic  constants  of  germanium  between  300  and  3°K  + 

James  A Dams.  Jr  and  M A.  Hrea/eale 
Department  of  Physics.  The  University  of  Tennessee.  Knoxville.  Tennessee  i 7916 
(Received  3 November  1975) 

The  ultrasonic  harmonic-generation  technique  has  been  improved  to  the  point  that  both  amplitude  and  phase 
of  the  second  harmonic  of  an  initially  sinusoidal  ultrasonic  wave  in  a solid  can  be  measured  between  room 
temperature  and  3’K  By  measuring  along  (he  principal  crystallographic  directions,  we  have  been  able  to 
determine  temperature  dependence  of  linear  combinations  of  third-order  clastic  (TOE)  constants  of 
germanium.  Between  room  temperature  and  77 *K  the  magnitude  of  the  TOE  constants  does  not  vary  greatly 
as  a function  of  temperature  Between  77  and  3*K.  C,(,  changes  by  +■  i%.  (C,,,  -f  4CtM)  changes  by  16%. 
and  (6C,„  + Cm  8Cllt)  changes  by  — 150%  All  of  these  combinations  of  TOE  constants  are  negative  at 
77*K.  but  below  7*K  the  combination  (6CI44  C..|  •+■  8C,H)  is  positive  Temperature  dependence  of 

additional  combinations  is  inferred  by  assuming  that  the  "Anderson-Ciruncisen  parameter"  is  independent  of 
temperature. 


i 


I INTRODUCTION 

The  present  work  ts  an  extension  to  very  low 
temperatures  ol  previously  established  techniques 
for  measuring  third-order  elastic  (TOE)  constants. 
The  technique  involves  the  measurement  of  the 
distortion  of  an  ultrasonic  wave  as  it  propagates 
through  the  solid.  This  distortion  was  observed 
in  polycrystalline  aluminum.  1 and  in  several  sin- 
gle crystals2  with  quartz  transducers.  The  mea- 
surement technique  was  refined  with  the  develop- 
ment of  a capacitive  receiver.1  This  refinement, 
and  subsequent  ones, 4 allowed  absolute  amplitude 
measurements  to  be  made  at  different  tempera- 
tures. Such  measurements  have  been  used  by  a 
number  of  investigators5  to  calculate  combinations 
of  truly  adiabatic  TOE  constants.  Yost  and  Brea- 
zeale*  combined  results  from  using  this  tech- 
nique with  those  of  Dunham  and  Huntington7  to  ob- 
tain the  first  complete  set  of  truly  adiabatic  TOE 
constants  (of  fused  silica)  at  room  temperature. 
TOE  constants  of  copper9  have  been  measured 
between  room  temperature  and  80° K.  More  re- 
cently we  reported®  data  on  germanium  between 
room  temperature  and  80  ’K.  (This  publication* 

Is  hereinafter  referred  to  as  I. ) The  present 
measurements  were  undertaken  to  extend  the  data 
on  TOE  constants  on  germanium  to  lower  tempera- 
tures in  an  attempi  to  answer  a fundamental  ques- 
tion about  the  relationship  between  the  clastic 
and  the  tcrmal  properties. 

The  thermal  expansion  of  germanium  is  anom- 
alous lielow  80  K. 10  A similar  anomaly  has  been 
observed  in  all  other  diamond  lattice  materials 
measured  to  date.11  At  temperatures  below  ap- 
proximately 0.  2 0n  where  is  the  Debye  tem- 
perature at  0 K,  the  thermal  expansion  changes 
from  positive  to  negative  and  then  becomes  posi- 
tive again  at  even  lower  temperatures.  Calcula- 
tions of  the  thermal  expansion  from  the  second- 


order  elastic  (SOE)  constants  and  TOE  con- 
stants12*11,14  have  not  predicted  the  measured 
anomalous  behavior  at  low  temperatures.  How- 
ever, these  calculations  were  based  on  the  assump- 
tion that  the  TOE  constants  are  not  dependent  on 
temperature  (in  lieu  of  actual  data  of  TOE  con- 
stant behavior  at  low  temperature).  We  chose 
germanium  for  these  low-temperature  measure- 
ments to  examine  the  validity  of  this  assumption. 

II  EXPERIMENTAL  TECHNIQUE 

Pure  mode  propagation  is  possible  for  a longi- 
tudinal ultrasonic  wave  in  the  three  principal 
directions  of  a cupic  crystal.  In  these  special 
directions,  the  wave  equation  takes  on  the  follow- 
ing form15: 

Port--  K,  ♦ (3/v,  . A-,)  (1) 

where  A2  and  A',  are  combinations  of  SOE  and  TOE 
constants,  respectively!  given  in  Table  I. 

If  there  ts  an  initially  sinusoidal  disturbance 
at  a - 0,  this  equation  has  the  solution 

n- Ay  sin (bn  - tot)  - [(3A'j  ♦ A's)  /8 Kt  | 

x A*b*a  cosMba  - tot),  (2) 

where  b is  the  propagation  constant  2v/\,  u is  the 


TABLE  I.  K,  and  K 3 lor  principal  directions  of  a 
cubic  crystal. 


Direction 

(100 1 

t'.l 

till 

(110 1 

• C,J  » 3 (.44) 

■i(c,„  . JC'ijj  • 1 

dill 

iff,,  • 2C'||  ♦ 4t'||) 

£«r,„  ♦ tit'iij  * 1 

12C  ||| 
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propagation  distance  In  the  sample,  and  A,  Is  the 
amplitude  of  the  fundamental. 

Thus  the  Initially  sinusoidal  wave  becomes  dis- 
torted: a second  harmonic  is  generated  and  throws 
linearly  with  propagation  distance  in  the  sample. 
The  amplitude  A2  of  this  second  harmonic  is  given 
by 

/t;  = - [<3A'2  » A',)/8A,  (3) 

This  harmonic  distortion  is  the  basis  of  the  non- 
linearity measurements  in  the  present  work.  The 
parameter  A',  may  be  calculated  by 

K2*  -3A8(/3*  1),  (4) 


0*HAt/A])(l/kzn),  (5) 

and  the  parameter  Kz  may  be  found  from  the  rela- 
tion 

Ki  * An'1*.  (6) 

where  r is  the  velocity  of  sound  In  the  appropriate 
direction  in  the  medium. 

A Rmim -lemperaiurr  mriturrmnili 

A block  diaitram  of  the  apparatus  for  making 
room -temperature  measurements  Is  shown  in 
Flu.  1.  A radio-frequency  pulse  of  approximately 
30  MHz  is  applied  to  an  x-cut  quartz  transducer 
which  Is  bonded  to  one  end  of  the  sample.  The 
ultrasonic  wave  which  propagates  through  the  sam- 
ple is  detected  at  the  other  end  by  a capacitive 
receiver.  In  addition,  provision  Is  made  to  place 
a substitutional  signal  on  the  capacitive  receiver 
so  that  accurate  amplitude  measurements  can  be 
made. 


The  air-gap  capacitive  receiver  used  is  essen- 
tially the  same  as  that  described  in  detail  by 
Gauster  and  Breazeale. J A 1.016-cm  diam  elec- 
trode is  placed  approximately  5 pm  from  the 
sample  (both  surfaces  being  optically  Hat/,  and 
a dc  bias  on  the  order  of  150  V is  applied  across 
the  gap  through  a large  resistor  (approximately 
1 MR). 

The  method  of  Introducing  the  substitutional 
signal  differs  somewhat  from  that  previously 
used  by  Peters,  Breazeale,  and  Pare*  and  in  I. 
Previously  the  capacitive  receiver  itself  had  been 
removed  from  the  circuit  during  the  calibration 
procedures,  and  the  calibrating  signal  was  intro- 
duced through  a substitutional  capacitor,  so  that 
the  equivalent  circuit  of  the  capacitive  receiver 
(an  ideal  voltage  generator  in  series  with  the  ca- 
pacitance of  the  receiver)  was  replaced  by  a real 
generator  (whose  voltage  can  be  measured  with  a 
voltmeter)  in  series  with  the  substitutional  ca- 
pacitor. (A  second  substitutional  capacitor  cor- 
rected for  stray  capacitance  from  the  receiver  to 
ground. ) 

In  the  present  measurements,  no  circuit  ele- 
ments were  substituted  during  the  calibration.  A 
substitutional  signal  which  gave  the  same  output 
as  the  ultrasonic  signal  was  introduced  across 
the  capacitive  receiver  in  such  a way  that  the 
current  of  the  substitutional  signal  could  be  mea- 
sured. This  current  i can  be  related  to  the  am- 
plitude A of  the  ultrasonic  wave  by 

i - 2A  V>u>c/a,  (7) 

where  V¥  is  the  dc  bias  on  the  receiver,  is  the 
angular  frequency,  r is  the  capacitance  of  the  re- 
ceiver. and  % is  the  spacing  of  the  capacitive  re- 
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cetver  plates. 

The  signal  from  the  capacitive  receiver  is 
taken  to  either  a 30-  or  60- MHz  bandpass  ampli- 
fier. This  amplified  signal  is  detected  and  taken 
to  a boxcar  integrator.  The  boxcar  Integrator 
selects  a portion  of  the  first  echo  and  gives  a sig- 
nal out  which  is  proportional  to  the  time  average 
of  that  portion  of  the  echo.  The  germanium  sam- 
ples. being  nonconductors,  do  not  sufficiently 
shield  the  rf  pulse  at  the  transducer  from  the  re- 
ceiver. so  that  a 30  MHz  pulse  is  received  at  the 
receiver  (by  radiation)  before  the  ultrasonic  wave 
reaches  the  receiver.  This  "feedthrough"  pulse 
can  be  much  larger  than  the  echo  of  interest  and. 
although  it  can  be  separated  in  time  from  the  mea- 
sured echo,  it  can  overload  the  30  MHz  amplifier 
and  the  boxcar  integrator.  This  problem  was 
solv»*d  by  gating  the  signal  from  the  capacitive 
receiver  so  that  only  the  echo  of  interest  was 
passed. 

Alter  the  fundamental  and  second-harmonic 
signals  have  been  measured  by  the  boxcar  inte- 
grator. the  continuous  wave  (cw>  substitutional 
signal  is  Introduced  at  the  capacitive  receiver. 

The  60  MHz  substitutional  signal  Is  derived  by 
doubling  the  30  MHz  signal  with  a ring  bridge 
mixer.  Both  the  30  and  GO  MHz  substitutional 
signals  are  filtered  to  ensure  spectra)  purity. 
These  two  signals  are  adjusted  to  give  the  same 
output  troni  the  boxcar  Integrator  as  the  ultra- 
sonic signals  and  are  then  measured  with  an  rf 
voltmeter.  From  these  measurements  and  a 
nowlerjgp  of  the  circuit  impedances  (which  were 
obtained  at  each  frequency  used  bv  measurements 


with  a vector  voltmeter' . the  current  of  the  sub- 
stitutional signal  can  be  calculated. 

B S«‘r»i>iJ  onJrr  riavlir  n>n\unl\ 

The  SOE  constants  were  determined  (rom  room- 
temperature  velocity  measurements  and  the  re- 
sults were  compared  with  the  data  of  McSklmin 
and  Adreatch16  and  McSkimln. 11  The  agreement 
Was  within  experimental  error.  The  extensive 
tabulated  data  of  McSklmin17  were  used  in  the 
calculations  of  the  cryogenic  data. 

C.  Cryogenic  3|jparalu\ 

The  electonic  apparatus  for  making  measure- 
ments at  low  temperature  is  essentially  the  same 
as  that  for  measuring  room  temperature. 

Since  relative  measurements  can  be  made,  it  is 
not  necessary  to  calibrate  at  each  temperature. 
Thus,  a reference  voltage  is  not  necessary. 

As  shown  in  Fig.  2,  the  sample  holder  is  en- 
closed in  a stainless-steel  can.  This  can  shown 
Is  surrounded  by  another  can.  The  space  between 
the  cans  is  evacuated  to  provide  an  insulating 
jac'cet  around  the  inner  can,  and  the  cans  are  pol- 
ished to  reduce  radiation  losses.  The  cans  are 
Supported  by  three  thin-walled  cupro-nickel  tubes. 
Two  of  the  tubes  have  a smaller  cupro-nickel 
tube  centered  Inside  to  provide  a 50- f?  coaxial 
transmission  line  and  are  also  used  as  vacuum 
lines.  The  other  tube  houses  the  temperature 
Sensor  and  heater  leads,  and  also  is  used  to  sup- 
ply a pressure  to  the  inner  can. 

The  capacitive  receiver  is  similar  to  that  in 
the  room  temperature  apparatus,  with  the  addi- 
tional feature  that  the  gap  spacing  can  be  con- 
trolled pneumatically. 4 The  lapped  ring  against 
which  the  sample  rests  Is  undercut  to  make  it  a 
flexible  dlnphram  approximately  0 038  cm  thick. 
The  space  above  the  sample  is  evacuated  and  the 
capacitor  gap  spacing  is  adjusted  by  regulating 
the  pressure  in  the  inner  can. 

The  entire  apparatus  is  suspended  Inside  a 
standard  helium  research  dewar.  The  dewar  can 
be  sealed  and  temperatures  below  the  room- 
pressure  boiling  point  of  the  coolant  may  be  ob- 
tained bv  pumping  on  the  coolant  (either  liquid 
helium  or  nitrogen). 

The  tomporaturo  in  the  Inner  can  Is  controlled 
by  an  electric  resistance  heater  connected  to  a 
commercial  temperature  controller,  and  can  be 
varied  continuously  from  approximately  .3  K to 
room  temperature. 

I)  Cryogenic  itniilinrariiy  measurement1* 

Only  relative  measurements  need  to  be  taken 
with  the  cryogenic  apparatus.  One  firsts  adjusts 
Hie  driving  signal  to  the  quartz  transducer  so  that 
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the  fundamental  ultrasonic  wave  received  at  the 
capacitive  receiver  is  the  same  at  each  tempera- 
ture. The  second  harmonic  is  then  measured  bv 
a slide- back  technique:  the  dc  bias  voltage  on  the 
capacitive  receiver  is  adjusted  so  that  the  elec- 
trical signal  coming  from  the  capacitive  receiver 
is  the  same  at  all  temperatures.  The  second- 
harmonic  amplitude  .U  of  the  ultrasonic  wave  at 
two  different  temperatures  (7*,  and  Tt)  are  then 
related  by 

At(Tt)  At(Tx)  V^Tx)/VtlTt).  (8) 

The  electrical  feedthrough  is  worse  in  the  cryo- 
genic apparatus  than  it  is  in  the  room  tempera- 
ture apparatus.  Therelore.  a 60- MHz  bandpass 
lllter  Is  used  between  the  capacitive  receiver  and 
the  amplifier  to  attenuate  the  30-MHz  feedthrough 
while  passing  the  60- MHz  signal. 

During  the  course  of  our  measurements  it  be- 
came necessary  to  determine  not  only  the  magni- 
tude of  the  second -harmonic  amplitude,  but  also 
its  sign.  This  resulted  from  the  fact  that  A,  were 
observed  to  change  enough  that  they  conceivably 
rould  have  gone  through  zero.  Such  a condition 
would  have  l>een  indicated  by  a phase  shift  of  180’ 
in  the  second  harmonic.  For  this  reason  we  per- 
fected a phas e- sensitive  detector  to  measure  the 
sign  of  the  second  harmonic.  (Phase- sensitive 
detectors  which  operate  at  60  MHz  are  not  com- 
mercially available  to  our  knowledge. ) 

A block  diagram  of  the  phase  sensitive  detector 
system  is  shown  in  Fig.  3.  A signal  from  the 
capacitive  receiver  Is  fed  through  a power  splitter 
to  a 30-MHz  amplifier  and  a GO- MHz  amplifier. 
The  output  of  the  30-MHz  amplifier  Is  fed  through 
a continuously  variable  phase  shifter  to  a fre- 
quency doubler.  This  frequency  doubled  30-MHz 
pulse  echo  train  and  the  second- harmonic  pulse 
echo  train  are  fed  into  the  phase- sensitive  detec- 
tor whose  output  is  proportional  to  the  cosine  ol 
•he  phase  angle  between  the  two  60- MHz  signals. 
The  phase  shifter  allows  the  phase  angle  between 
the  two  signals  to  be  set  to  zero  initially.  If  the 
sign  of  the  second- harmonic  amplitude  changes 
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as  the  temperature  is  changed,  the  output  of  the 
phase-sensitive  detector  is  Inverted, 

Using  this  phase  sensitive  detector  system,  we 
were  able  to  make  an  unequivocal  assignment  of 
the  phase  to  Hie  second- harmonic  amplitude  used 
in  the  calculations  of  the  linear  combinations  of 
TOE  constants  we  call  A'*. 

I.  Samples 

The  samples  used  In  this  work,  were  the  same 
ones  used  in  1.  These  single  crystals  were  3.7803. 
4.8242.  and  3.4803  cm  in  the  | i 00 1,  illOj.  and 
fill | orientations,  respectively.  (The  lengths  are 
slightly  cnanged  since  used  in  I because*  ot  lapping 
and  recoating  of  the  sample  ends. ) The  sample 
ends  were  lapped  and  polished  until  optically  flat 
and  parallel  to  within  12  sec  of  arc.  The  ends 
were  then  made  conductive  by  evaporating  a cop- 
per coating  approximation  1000  A thick  onto  them. 

Ill  RFSLLTS 

A Room  temperature  nonlinearity  measurement*  lor  Ge 

The  absolute  amplitude  measurements  at  room 
temperature  are  given  in  Fig.  4.  The  straight 
lines  in  the  figures  are  least-squares  fits  to  the 
data  points.  It  can  be  seen  that  the  data  will  fit 
very  well  with  a straight  line,  however,  this  line 
does  not  pass  through  the  origin  because  of  residu- 
al noise.  The  slope  of  the  line  is  (In  the  least- 
squares  sense)  the  best  fit  to  the  value  of  the  quan- 


moiuc  amplitude  plotted  an  a (unction  ot  the  M|uurc  ol  (ho 
fundamental  amplitude. 
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TABLE  II.  Measured  room-temperature  values  <>t  p and  A.. 
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TEMPERATURE  f*K) 

FIG.  6.  Temperature  dependence  of  the  combinations 
of  TOE  constants  A,  for  the  three  principal  directions  in 
Germanium. 

of  temperature.  The  first  data  runs,  shown  by 
the  open  symbols  were  taken  without  use  of  the 
phase- sensitive  detector.  The  curves  for  the 
f 100 J and  f 111  J directions  are  considered  to  be 
too  smooth  to  admit  a negative  phase  of  the  sec- 
ond harmonic  from  which  the  data  were  calcu- 
lated. The  data  in  the  fll0|  direction  were  not 


TEMPERATURE  (°K) 

FIG.  7.  remperatwre  «k*|«enik,iK*e  of  the  simplest  com- 
Mfialiom  if  TOE  i onstanl*  -Ahleh  can  «ab  ulatcd  front 
the  h ol  (.<  i minium. 


100  200  300 

TEMPEPATORe  (°K  ) 

FIG.  1.  Temperature  dependence  of  certain  TOE 


constants  combinations  calculated  from  our  data  under 
the  assumption  that  the  "Anderson-Griineisen  param- 
eter" and  Cia  are  temperature  independent. 


as  smooth.  The  initial  data  between  approximately 
15  and  50  JK  had  shown  a peak  or  dip  depending 
on  the  interpretation  of  the  phase.  We  also  noted 
that  a thermal  gradient  might  have  been  respon- 
sible for  this  behavior.  The  apparatus  was  modi- 
fied to  eliminate  thermal  gradients  and  the  phase- 
sensvUve  detector  was  added.  The  data  shown  as 
solid  circles  were  taken  with  the  modified  appara- 
tus. There  is  no  ambiguity  about  the  sign  of  the 
second-harmonic  amplitude  from  which  these  data 
were  calculated.  The  agreement  between  the 
two  sets  of  data  In  the  low-temperature  range 
reassured  us  about  the  correctness  of  these  data 
and  allowed  us  to  eliminate  erroneous  data  points 
from  this  plot. 

Down  to  liquid- nitrogen  temperature,  these 
data  reproduce  the  results  of  I to  well  within  TTc. 
The  data  for  the  (l00|  direction  can  be  fit  very 
well  with  a straight  line  down  to  3r  K.  The  data 
for  the  f 111  | direction  is  essentially  a straight 
line  with  a very  slight  upturn  at  very  low  tem- 
peratures. The  data  for  the  1 1 1 0 i direction  are 
the  only  ones  which  show  a marked  deviation  from 
a straight  line.  At  low  temperatures  there  is  a 
definite  downturn  of  the  data. 

These  linear  combinations  of  the  TOE  constants 
are  not  the  simplest  possible  combinations  avail- 
able from  our  data.  For  example.  Cm  Is  con- 
tained in  each  of  these  curves,  and  therefore  can 
be  subtracted  out.  ''roceedlng  in  this  fashion,  we 
have  evaluated  the  simplest  combinations  ot  TOE 
constants  shown  in  Fig.  7. 

As  can  Ik*  seen,  there  is  not  a great  variation 
of  these  combinations  ol  third-order  elastic  con- 
stants down  to  77  K.  But  Iwtween  77  and  3 K 
there  is  a considerable  variation  in  some  of  them. 
In  particular,  the  combination  It’,,,  ♦ 6C1M  ♦ SC,*) 
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FIG.  a.  Temperature  dependence  of  certain  TOE  con- 
stant combination.-,  calculated  trom  our  data  assuming 
that  lx>th  the  "Anderson-GrUneisen  parameter  " and  L\a 
are  ienq»erature  independent. 


crosses  zero  at  7 K and  becomes  positive  below 
that  temperature.  Since  these  are  the  TOE  con- 
stants which  would  play  a significant  role  in  the 
transverse  modes,  we  leel  that  the  lack  of  agree- 
ment between  theoretical  calculation  of  thermal 
expansion  from  elastic  data12*13*14  and  the  mea- 
sured thermal  expansion  might  be  attributed  to 
the  theoretical  assumption  that  the  TOE  constants 
are  independent  of  temperature. 

Ill  ISOLATION  OF  ADDITIONAL  TOE  CONSTANTS 

Further  reduction  of  the  data  can  be  made  if  one 
makes  the  assumption  that  the  " Anderson-Griin- 
elsen  parameter"1*  is  not  a function  of  temperature.  19 
This  assumption  may  be  more  valid  than  the  as- 
sumption that  the  TOE  constants  are  independent 
of  temperature.  Therefore,  using  it  to  isolate 
additional  TOE  constants  from  our  data  may  prove 
to  be  instructive. 

Using  the  data  from  our  experiment,  the  data 


of  McSkimm  and  Andreatch14  [combination  of  TOE 
constants  at  room  temperature),  and  the  data  of 
MeSklmin17  for  the  second-order  constants  Cu 
ar.d  C 12  we  can  calculate  the  “Anderson-Gruneisen  ’ 
parameter  by  Rao’s  formulation20: 

^ ; ~ 1 ~ (Cm  ♦ 6C  u2  ♦ 2C12j) /3(C”u  ♦ 2CI;I.  ( 10) 

If  5 is  indeed  constant  as  a function  of  tempera- 
ture. then  the  temperature  dependence  of  the 
combination  3C,,2  + C123  can  be  obtained  from  our 
data.  Under  this  assumption,  this  combination 
of  TOE  constants  shows  virtually  no  temperature 
dependence  (see  Fig.  8).  If  it  is  further  assumed 
that  C,23  is  not  a function  of  temperature  (and 
equal  to  the  room  temperature  value  obtained  from 
the  data  of  our  experiment  and  McSkimin  and 
Andreatch),  additional  TOE  constants  mav  be  cal- 
culated. These  calculated  constants  are  given 
in  Fig.  9 and  are  designated  by  enclosing  the  con- 
stants in  quotation  marks  to  specify  that  these 
calculations  are  based  on  the  two  assumptions 
given  above. 

In  conclusion,  the  measurements  indicate  that 
there  is  a definite  temperature  dependence  of  some 
of  the  TOE  constants.  Since  we  cannot  isolate 
all  of  the  TOE  constants,  it  is  not  feasible  at 
present  to  determine  the  elfect  of  this  tempera- 
ture dependence  on  the  Griineisen  parameter  cal- 
culation of  Brugger  and  Fritz,  14  However,  our 
assumption  that  the  Anderson-Gruneisen  param- 
eter is  independent  of  temperature  has  allowed  us 
to  plot  most  of  the  TOE  constants  as  a function  of 
temperature.  The  validity  of  this  assumption, 
and  hence  of  the  curves  drawn  to  Figs.  8 and  9. 
is  open  to  question;  however,  one  can  observe 
that,  to  the  extent  these  curves  are  valid,  the 
greatest  temperature  dependence  seems  to  be 
found  in  C1M  and  3Cl44  ♦ 4 C456. 
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ULTRASONIC  NONLINEARITY  PARAMETERS  OF  FUSED  SILICA  BETWEEN 
3SK  and  300SK 
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Cantrell,  J.  H.  The  University  of  Tennessee 

Knoxville,  Tennessee  37916  USA 


The  nonlinearity  parameters  of  fused  silica  have  been  measured  between  room 
temperature  and  3*  Kelvin  by  the  harmonic  generation  technique  (1] . 

A 30  MHz  ultrasonic  pulse  is  generated  in  samples  of  fused  silica  which  differ 
in  hydroxyl  ion  content  and  in  homogeneity.  As  the  ultrasonic  wave  propagates 
through  the  sample  it  generates  harmonics  of  the  30  MHz  fundamental.  By  use  of  a 
capacitive  receiving  transducer,  amplitudes  of  both  the  fundamental  and  the  second 
harmonic  are  measured.  Typical  values  of  the  fundamental  and  second  harmonic 
amplitudes  are,  respectively,  2.2  Angstroms  and  9.3  x 10"3  Angstroms.  These  ampli- 
tudes are  used  to  calculate  the  nonlinearity  parameter  (2) 
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and  hence  the  third-order  elastic  constant  C....  Room  temperature  values  of  these 
quantities  for  four  different  samples  are  given  in  Table  1. 

As  the  temperature  is  lowered  to  3 *K,  Cj, j changes  by  16%  for  Suprasil  Kl,  and 
less  than  this  for  the  other  samples.  The  variation  of  B is  7%  for  Suprasil  2,  and 
less  than  this  foT  the  other  samples.  Since  0 can  be  related  to  the  longitudinal 
mode  strain  Grtlneisen  parameter  Yjj,  this  means  that  temperature  dependence  of  the 
total  Grtlneisen  parameter  as  great  (3)  as  800%  must  come  from  sources  other  than 
variations  of  the  third-order  elastic  constant  C....  [Research  sponsored  by  the 
Office  of  Naval  Research.] 

TABLE  I 


Sample  B Cn  (1011  dynes/cm2)  Cm  (1012  dynes/cm2) 


Suprasil  Ml 

-3.86  ± 

0.071 

7.8059  ± 

0.0029 

6.70 

1 0.124 

Suprasil  W2 

-3.93  t 

0.034 

7.8078  ± 

0.0020 

6.85 

1 0.050 

Suprasil  1 

-4.14  ± 

0.041 

7.7633  t 

0.0016 

7.31 

t 0.073 

Suprasil  2 

-3.90  ± 

0.039 

7.7633  ± 

0.0014 

6.75 

± 0.067 
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Ultrasonic  investigation  of  the  nonlinearity  of  fused  silica  for  different  hydroxyl-ion  contents 
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John  H.  Cantrell,  Jr.*  and  M.  A.  Breazeale 
Department  of  Physics,  The  University  of  Tennessee.  Knoxville,  Tennessee  37916 
(Received  19  September  1977) 

The  harmonic-generation  technique  was  used  to  measure  the  third-order  elastic  constants  Cm  of  four  types 
of  fused  silica.  The  samples  differed  in  amount  of  OH  content  and  directional  homogeneity.  The  results 
indicate  that  a relatively  large  OH  content  in  fused  silica  may  give  rise  to  a relatively  greater  value  of  C,n 
and  that  CU1  is  a weak  function  of  temperature  regardless  of  OH  content.  The  relationship  between  the 
experimentally  defined  ultrasonic  nonlinearity  parameter  0=  — (3  C, , + Cm)/3CM  and  the  Griineisen 
parameter  y also  has  been  determined.  The  measured  temperature  dependence  of  0 indicates  that  in  the 
Debye  model  elastic  constants  other  than  Cn  and  C|M  must  dominate  the  strong  temperature  dependence  of 
the  Griineisen  parameter  for  fused  silica  at  low  temperatures. 


I.  INTRODUCTION 

A number  of  the  physical  properties  of  solids 
result  from  the  fact  that  solids  are  inherently  non- 
linear. Among  these  properties  are  thermal  ex- 
pansion, attenuation  of  high-frequency  sound 
waves,  heat  conduction,  and  wave-form  distortion 
of  ultrasonic  waves  passing  through  a solid.  These 
nonlinear  effects  arise  because  of  anharmonicity 
of  the  interatomic  potential. 

The  earliest  investigations  of  anharmonicity  in 
solids  were  pursued  from  more  purely  thermody- 
namic experimentation.  For  example,  Fizeau1 
was  among  the  earliest  investigators  to  measure 
thermal  expansion. 

Gr&neisen5  developed  the  relationship  between 
the  thermal-expansion  coefficient  and  a parameter 
which  was  assumed  to  be  temperature  independent 
(the  Griineisen  parameter  y).  Barron5'4  and 
Sheard5  were  among  the  first  to  make  calculations 
relating  the  Griineisen  parameter  to  elastic  data. 
They  limited  their  calculations  to  the  cases  of  high 
and  low  temperatures.  Collins,*  Schuele  and 
Smith,7  Brugge r and  Fritz,*  and  Gerlich*  have 
made  more  extensive  calculations  to  determine  the 
temperature  dependence  of  the  Ouneisen  parame- 
ter from  elastic  data  using  the  quasiharmonic 
Debye  model  of  solids.  Thus,  the  Ouneisen  con- 
stant is  a convenient  link  between  thermal-expan- 
sion data  and  elastic  data. 

The  agreement  of  the  calculated  dependence  of 
the  Oiinelsen  .iarameter  as  a function  of  tempera- 
ture with  values  measured  from  thermal  expansion 
has  been  limited  because  of  the  lack  of  availability 
of  elastic  data  [particularly  third-order  elastic 
(TOE)  constants]  as  functions  of  temperature.  For 
many  materials  such  as  Cu,  Ag,  Au,  Al,  and  Na 
qualitative  agreement  exists,  but  some  materials 
exhibit  an  anomalous  behavior  of  the  thermal -ex- 
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pansion  coefficient  as  a function  of  temperature. 
One  such  material  is  fused  silica.  Fused  silica 
has  a positive  thermal-expansion  coefficient  down 
to  a temperature  of  approximately  0.3  of  its  Debye 
temperature.  There  the  thermal-expansion  coef- 
ficient becomes  negative  and  continues  to  grow 
more  negative  as  the  temperature  is  lowered.10 
This  gives  rise  to  the  anomalous  behavior  of  the 
Griineisen  parameter  y shown  in  Fig.  1,  where  it 
is  seen  that  y also  becomes  negative  around  0.3  of 
the  Debye  temperature  and  continues  to  become 
more  negative  as  the  temperature  approaches  zero. 
This  is  in  contrast  with  germanium,  for  which  y 
levels  off  as  the  temperature  approaches  zero  and 
with  quartz  which  maintains  a positive  y for  all 
temperatures.  Several  investigations5-®' 11-15  indi- 
cate that  at  least  part  of  the  explanation  for  the  be- 
havior of  fused  silica  may  be  due  to  the  dominance 
of  transverse-acoustical  modes  at  these  tempera- 
tures; however,  this  hypothesis  can  be  tested  only 
as  information  on  the  behavior  of  TOE  constants 
as  functions  of  temperature  is  collected.  Since  the 
TOE  constants  determine  how  a large  amplitude 
ultrasonic  wave  distorts  as  it  propagates  through 
the  solid,  TOE  constants  can  be  calculated  from 
ultrasonic  wave  distortion.  We  have  chosen  to 
measure  this  wave-form  distortion  in  order  to 
understand  more  fully  the  unusual  behavior  of 
fused  silica. 

Ultrasonic  wave  distortion  was  initially  observed 
in  polycrystalline  aluminum1*  and  in  several  single 
crystals17  with  quartz  transducers.  Absolute  am- 
plitude measurements  became  possible  with  the 
development  of  the  capacitive  receiver,1®  and  mea- 
surements at  different  temperatures  became  possi- 
ble with  a later  refinement.1®  Several  investiga- 
tors5045 have  used  this  method  to  calculate  com- 
binations of  truly  adiabatic  TOE  constants.  Yost 
and  Breazeale55  combined  the  results  of  this  meth- 
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FIG.  1.  Temperature  dependence  of  the  Griinelsen  parameter  for  quartz,  germanium,  and  spectrosll  fused 
silica  [from  G.  K.  White,  Cryogenic  4,  2 (1964)1 . 

od  with  those  of  Dunham  and  Huntington24  to  obtain 
the  first  complete  set  of  purely  adiabatic  TOE  con- 
stants (of  fused  silica)  at  room  temperature.  Com 
binations  of  TOE  constants  of  copper25  and  german 
ium26  have  been  measured  between  room  tempera- 
ture and  77  °K.  More  recently,  the  measurements 
on  germanium  were  extended  to  lower  tempera- 
tures (3  TO.27 

The  purpose  of  the  present  investigation  is  to  ex- 
plore the  results  of  anharmonicity  more  fully.  To 
do  this,  we  have  measured  the  second  harmonics 
of  an  initially  sinusoidal  ultrasonic  wave  in  four 
different  samples  of  fused  silica  between  room 
temperature  and  3 °K.  From  these  data  we  have 
calculated  the  adiabatic  TOE  constant  Cm  as  a 
function  of  temperature.  Also,  we  have  plotted  the 
experimentally  defined  nonlinearity  parameter 
0 = -(3CU  +Clu/3CU),  because,  as  we  show,  there 
is  a relationship  between  0 and  a Griineisen  pa- 
rameter. The  implications  of  the  present  results 
in  light  of  the  temperature  dependence  of  the 
Griineisen  parameter  are  discussed. 

n.  THEORY 

A.  Ultrasonic  nonlinearity  in  isotropic  solids 

For  an  isotropic  solid  pure  mode  propagation  is 
possible  for  a longitudinal  ultrasonic  wave  in  any 


direction.  In  this  special  case  the  nonlinear  wave 
equation  takes  the  form28 

*«-cu£i*<3c, a, 

where  p0  is  the  mass  density  of  the  undisturbed 
medium,  u is  the  particle  displacement,  a is  the 
Lagrangian  (laboratory)  coordinate  in  any  direc- 
tion, and  C„  and  CU1  are  second-  and  third-order 
elastic  constants,  respectively. 

Assuming  an  initially  pure  sinusoidal  disturbance 
at  a =0,  the  solution  to  Eq.  (1)  may  be  written2* 

u=A,sin(*a  -u if)  -[(3C„  +Cm)/8Cll]A}*2 

xacos2(fea-wf),  (2) 

where  k is  the  propagation  constant  2ir/A,  and  u> 
is  the  angular  frequency.  Hence,  an  initial  sinus- 
oidal disturbance  of  amplitude  A,  distorts  linearly 
with  propagation  distance  a and  generates  a second 
harmonic  of  amplitude 

A1  = 4(3ClI+C1I1)/8C11]A;*2a.  (3) 

Solving  Eq.  (3)  for  Clu  yields 

Cm  = -3CU[ } W,/A2)(  1/k  2a)  + 1] . (4) 

All  quantities  on  the  right-hand  side  of  this  equa- 
tion can  be  measured  and  C...  can  be  determined. 
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The  quantity 

l(A,/A*l)(\/k*a)=0  (5) 

in  Eq.  (4)  is  a significant  quantity  in  this  study.  If 
0 is  zero,  then  A2  is  zero  and  no  second  har- 
monic is  generated.  Since  the  generation  of  a 
second  harmonic  is  a direct  measure  of  the  non- 
linearity of  a solid  then  (3  may  be  referred  to  as 
the  ultrasonic  nonlinearity  parameter.  For  iso- 
tropic solids  0 is  expressed  in  terms  of  the  elastic 
constants  as 


ffD  = (6nl/V0)l/,l  where  V0  is  the  volume  per  atom]. 
The  specific  heat  of  mode  ( p,N ) in  Eq.  (9)  is  given 
by 


C(p,ft)  = 


i*\e(p,N)/T]2 

{exp|e(/>,  A^)/7']  - l}2 


(10) 


where  0(p,N)  is  the  characteristic  Debye  tempera- 
ture of  mode  ( p,N ) given  by 

6(P,ft)  = (hqD/K)S(p,N),  (11) 


0 = -(3C ,,  ♦C,11)/3Cll . 


(6) 


H (iriineurn  parameter*  and  their  relationship  to  0 
I (iriinelsen  parameters 

The  anharmonicity  of  solids  also  can  be  investi- 
gated from  thermodynamic  measurements  (e.g., 
volume  expansivity,  variations  in  isothermal  com- 
pressibility as  a function  of  temperature,  etc.). 
Among  those  investigating  anharmonicity  from  these 
techniques  it  has  become  common  to  express  the 
results  in  terms  of  the  Griineisen  parameter  y, 
defined  by2 


y = a/KrCv  = a/KsC  p,  (7) 

where  a is  the  thermal  volume  expansivity,  KT 
and  Ks  are  the  isothermal  and  isentropic  compres- 
sibilities, respectively,  and  Cv  and  Cp  are  the 
isochoric  and  isobaric  heat  capacities,  respective- 
ly- 

A relationship  between  y and  the  elastic  constants 
can  be  found  from  calculations  based  on  the  as- 
sumption that8’7"8  (i)  the  material  behaves  like 
continuous  anisotropic  medium,  (ii)  the  nondisper- 
sive  Debye  model  of  specific  heats  is  valid,  and 
(ill)  the  generalized  Griineisen  parameters  depend 
on  temperature  only  through  the  lattice  dimensions. 

If  one  defines  the  volume  generalized  Griineisen 
parameter  for  the  mode  (/>,q)  by 


(8) 


where  p is  the  polarization  index  (p  - 1, 2, 3)  cor- 
responding to  the  three  acoustical  modes,  q is  the 
wave  vector,  V is  the  volume,  T is  the  tempera- 
ture, and  u>(/>,  q)  is  the  angular  frequency  of  the 
mode.  The  Griineisen  parameter  y may  then  be 
expressed  as8 


y = £ fdUy(p,N)C(p,N )j ^fdQC(p,N) , 


(9) 


where  the  integration  is  performed  over  the  spatial 
direction  N = q/|q|  in  the  irreducible  part  of  the 
Brillouin  zone  [in  our  case  a sphere  of  radius 


where  K is  the  Boltzmann  constant  and 

S(p,jJ)=(?^5))  (12) 

\ 8</  / ro 


is  the  elastic  wave  speed  of  that  mode. 

It  is  convenient  to  introduce  the  strain  general- 
ized Griineisen  parameters  ylk(p,  q)  defined  by 


5>  = "[<*>(/>,  q) 


/ 8u)(/>,q)\  ‘ 

\ Sjjy,  ' T\V0> 


(13) 


where  q is  the  Lagrangian  strain  tensor  and  the 
indices  j,k  = 1,  2,  3.  The  volume  generalized 
Griineisen  parameters  can  be  expressed  in  terms 
of  the  strain  generalized  Griineisen  parameters  by 


y(/>.q)=^-  T (14) 

where  KT  is  the  isothermal  compressibility  and 
S*jrt  are  the  isothermal  elastic  compliance  coeffi- 
cients. 

It  has  been  shown28  that  ylt(p,  q)  can  be  expressed 
in  terms  of  the  elastic  constants  by  (Einstein  nota- 
tion) 

rlll(p,N)=ll/2w(p,N)][2w(p,N)ulull 

(15) 

where 

«>(P,N)°CMmSmN„U'V0.  (16) 

ft  gives  the  direction  of  wave  propagation,  0 is  a 
unit  vector  along  the  direction  of  polarization  ap- 
propriate to  p,  and  the  subscripted  C’s  are  the 
second-  and  third-order  elastic  constants  in  the 
equations. 


2 Relationship  to  the  ultrasonic  nonlinearity  parameter 

For  pure  longitudinal  acoustic  modes  in  an  iso- 
tropic solid,  Eq.  (15)  reduces  to  the  set  of  rela- 
tions (using  Voigt50  notation) 

,yll  = -(3Cll+C111)/2Cll,  (17) 

iVjj  = iVSj  = -(C|i +CU1)/2CU  , (18) 

and 
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,ru  = 0,  i*J,  (19) 

where  the  prescript  / indicates  longitudinal-acous- 
tical modes.  From  Eq.  (6)  and  (17)  it  is  seen  that 
,y,,  is  related  directly  to  the  ultrasonic  nonlinear- 
ity parameter  0 by 

,ru  = i0.  (20) 

Hence,  one  of  the  strain  generalized  Griineisen 
parameters  is  measured  directly  by  the  present 
harmonic-generation  technique. 

For  an  isotropic  solid  Eq.  (9)  reduces  to  the 
form 

y = (y(C,  + 2y,C,)/(C1  + 2C,),  (21) 

where 

y,  = +i/3-(C11  + Clll)/3Cu  (22) 

is  the  longitudinal -mode  volume  generalized 
Griineisen  parameter  obtained  from  Eq.  (8),  and 
C,  and  C,  are  the  modal  specific  heats  of  Eq.  (10). 
y,  in  Eq.  (21)  is  the  transverse-mode  volume  gen- 
eralized Griineisen  parameter  and  is  not  easily 
expressed  in  terms  of  the  ultrasonic  nonlinearity 
parameter. 

The  significance  of  the  relationship  of  the  non- 
linearity parameter  0 to  the  Griineisen  parameter 
y is  now  established  for  the  case  of  isotropic  solids, 
ft  allows  one  to  ascertain  directly  the  contribution 
of  the  “longitudinal-mode”  elastic  constants  Cn 
andCIuto  the  temperature  dependence  of  the 
GHineisen  parameter  when  0 is  measured  as  a 
function  of  temperature. 

Ill  EXPERIMENTAL  TECHNIQUE 
A.  Apparatus  and  measurement  techniques 

The  apparatus  and  measurement  techniques  for 
determining  the  amplitudes  of  fundamental  and 
second -harmonic  waves  are  described  in  detail 
elsewhere.21  The  amplitudes  of  the  fundamental 
and  the  second-harmonic  waves  were  measured 
absolutely  in  each  sample  at  room  temperature. 
These  data  were  used  in  Eqs.  (4)  and  (5)  to  deter- 
mine room -temperature  values  of  Clu  and  0.  At 
all  other  temperatures  the  amplitudes  were  mea- 
sured relative  to  the  room -temperature  values, 
and  relative  values  of  Cln  and  0 were  determined. 


B Second -order  elastic  constants 

The  second-order  elastic  (SOE)  constant  C,,  ap- 
pearing in  Eq.  (4)  was  determined  from  room-tem- 
perature longitudinal  wave-velocity  measurements, 
since 

C„=p0^,  (23) 

where  p0  is  the  mass  density  of  the  unstrained 
solid.  The  measurements  were  made  with  capaci- 
tive transducers31  which  eliminated  the  problem 
of  bond  corrections.  The  temperature  dependence 
of  Cn  was  calculated  from  the  data  of  A.  Zarem- 
bowitch32  who  measured  the  SOE  constants  of  Pur- 
posil  fused  silica  as  a function  of  temperature. 

The  SOE  constants  of  the  Suprasil  fused-silica 
samples  used  in  these  experiments  are  assumed 
to  have  the  same  relative  temperature  dependence 
as  that  of  Purposil. 

C.  Samples 

The  fused-silica  samples  used  in  this  work  were 
manufactured  under  the  commercial  designation 
Suprasil.  All  samples  were  cylindrical,  approxi- 
mately 2.54  cm  in  diameter,  and  were  coated  on 
each  end  with  copper  to  a thickness  of  approxi- 
mately 1000  A.  The  length  of  each  sample,  impur- 
ity content,  and  the  measured  ultrasonic  wave  vel- 
ocities at  room  temperature  (27  °C)  are  given  in 
Table  I.  The  symbol  W indicates  that  the  sample 
has  only  5 ppm  of  OH,  while  the  unlettered  sam- 
ples have  OH  contents  of  1200  ppm.  The  numerical 
designation  1 or  2 indicates  the  degree  of  homo- 
geneity of  the  sample.  The  designation  1 indicates 
a guaranteed  strong  homogeneity  in  all  directions, 
whereas  the  designation  2 implies  a guaranteed 
strong  homogeneity  only  in  the  direction  perpen- 
dicular to  the  sample  surfaces. 

IV.  RESULTS 

A.  Room-temperature  measurements 

Accurate  room -temperature  values  of  the  non- 
linearity parameter  0 and  of  C1U  were  determined 
by  an  extrapolation  technique  which  minimized  the 
effect  of  attenuation  on  the  results.”  The  nonlin- 
earity parameter  0 was  measured  for  different 


TABLE  I.  Samples  used  in  this  study  and  pertinent  properties. 


Impurity  content  (ppm)  Sample  length  Longitudinal  ultrasonic 
Sample  OH  Fluorine  Chlorine  (cm)  velocity  (10*  cm/sec) 


Suprasil  W1 

5 

ma 

260 

1.2590  * 0.0003 

5.9526  *0.0016 

Suprasil  W2 

5 

uuM 

260 

1.2587  * 0.0002 

5.9533  *0.0011 

Suprasil  1 

1200 

130 

1.2593  *0.0002 

5.9363  * 0.0012 

Suprasil  2 

1200 

■n 

130 

1.2586  *0.0002 

5.9363  * 0.0011 
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FIG.  2.  Room-temperature  plots  of  0 as  a function  of 
the  fundamental  amplitude. 


fundamental  wave  amplitudes  A,  for  each  sample. 
The  resulting  values  of  0 were  plotted  as  a function 
of  .4,,  as  shown  in  Fig.  2.  The  solid  curves  are 
drawn  through  the  data  points  in  such  a way  that 
they  approach  A,  sO  with  a horizontal  tangent. 
Previous  analysis”  shows  that  this  extrapolation 
corrects  for  the  effect  of  attenuation.  As  can  be 
seen  the  maximum  magnitude  of  this  correction  is 
of  the  order  of  5%. 

The  results  of  the  extrapolated  room-tempera- 
ture measurements  of  0,  C,,,  and  Cm  are  given  in 
Table  II.  An  examination  of  the  table  reveals  that 
the  SOE  constants  of  the  fused-silica  samples  with 
high  OH  content  (Suprasil  1 and  Suprasil  2)  are 
lower  than  the  samples  with  low  OH  content 
(Suprasil  Wl  and  Suprasil  W 2).  These  results  are 
consistent  with  those  of  Hetherington  and  Jack” 
who  noted  also  that  an  increase  in  OH  content 


AND  M . A . II  K F.  A 7.  K.  A I.  K, 

caused  an  increase  in  mass  density  fluctuation  as 
much  as  3 parts  in  103  in  their  experiments. 
Primak3''  suggested  that  in  similar  experiments 
by  Fraser35  in  the  ultrasonic  range  of  frequencies 
the  results  reflect  the  effect  of  network  cleavages 
induced  by  OH  impurities. 

Table  II  shows  that  0 is  negative  for  all  fused- 
silica  samples.  This  means  that  A2  is  also  nega- 
tive and  that  CU1  is  positive.  Bains  and  Breazeale36 
have  shown,  using  a phase-sensitive  detector,  that 
the  harmonics  of  fused  silica  are  generated  out  of 
phase  with  those  generated  by  the  same  ultrasonic 
wave  in  copper.  They  conclude  that  the  positive 
C1M  of  fused  silica  means  that  the  solid  becomes 
softer  upon  compression  in  contrast  to  copper 
which  becomes  stiffer  upon  compression. 

According  to  Table  II,  Suprasil  1 has  a substan- 
tially higher  value  of  Clu  than  any  of  the  other 
samples  used  in  these  experiments.  It  is  tempting 
to  infer  that  this  is  because  of  the  high  OH  content 
of  the  sample  (1200  ppm  compared  to  5 ppm  for 
Suprasil  Wl  and  Suprasil  VV’2),  but  Suprasil  2 also 
has  a high  OH  content  (1200  ppm)  and  gives  a value 
of  CU1  which  lies  between  those  of  Suprasil  Wl  and 
Suprasil  W 2.  The  reason  for  this  is  not  completely 
understood  at  this  time.  A check  of  all  the  samples 
with  crossed  Polaroid  sheets  revealed  that  residual 
stresses  were  undetectable  by  this  technique. 
Hence,  residual  stresses  in  a particular  sample 
probably  are  not  the  cause  of  the  difference. 

A possible  explanation  may  lie  in  the  fact  that 
even  though  Suprasil  2 is  strongly  homogeneous 
in  the  direction  perpendicular  to  the  sample  faces, 
it  may  have  inhomogeneities  in  other  directions. 

If  mass  density  variations  existed  parallel  to  the 
sample  faces,  then  there  would  be  an  associated 
phase  variation  across  the  ultrasonic  wave  front. 
Such  a situation  leads  to  phase  heterodyning  which 
can  cause  significant  changes  in  the  amplitude  of 
the  electrical  signal  response.37  Since  this  hetero- 
dyning is  frequency  dependent,  fundamental  and 
second-harmonic  waves  would  be  affected  differ- 
ently. The  result  could  be  a decrease  in  the  mea- 
sured value  of  Cm  for  Suprasil  2.  If  this  were 
true,  then  one  could  consistently  associate  high 
C1H  values  with  high  OH  content.  Thus,  the  net- 
work cleavages  introduced  by  the  OH  impurities 


TABLE  n.  Measured  room-temperature  values  of  0,  Clt,  and  Ctll.* 

Sample  0 C,i  (1011  dyn/cm2)  C(II  (1012  dyn/cm2) 


Suprasil  Wl 
Suprasil  Wt 
Suprasil  1 
Suprasil  2 


-3.86  •0.071 
-3.93  i 0.034 
-4.14  i 0.041 
-3.90  • 0.039 


7.8059  i 0.0029 
7.8078  i 0.0020 
7.7633  • 0.0016 
7.7633  t 0.0014 


6.70  i 0.124 
6.85  • 0.058 
7.31  * 0.073 
6.75  • 0.067 


The  errors  listed  are  the  calculated  standard  deviations  (random  errors). 
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FIG.  3.  Temperature  dependence  of  0 for  the  fused- 
stllca  samples. 
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which  may  be  the  cause  of  lower  C,,  values  could 
also  be  the  cause  of  higher  C,,,  values. 

B.  Cryogenic  measurements  of  0 and  its  implications 
to  the  calculation  of  Griineisen’s  y 

In  making  relative  measurements  with  the  cryo- 
genic system,  the  quantity  0 is  determined  direct- 
ly27 | see  Eq.  (5)].  The  resulting  plots  of  0 versus 
temperature  in  Fig.  3 indicate  that  0 is  virtually 
temperature  independent.  This  means  that  the  ,yn 
component  of  the  strain  generalized  Griineisen  pa- 
rameter also  is  virtually  temperature  independent 
since  ,y„  = i0. 

Therefore,  the  temperature  variation  of  the  ex- 
perimentally determined  Griineisen  parameter  y 
for  fused  silica  shown  in  Fig.  1 must  come  from 
variations  in  the  volume  generalized  Griineisen 
parameter  [see  Eqs.  (21)  and  (22)]. 

It  is  apparent  from  Eq.  (21)  that  in  order  to  ob- 
tain a negative  value  of  the  Griineisen  parameter 
at  some  temperature  for  fused  silica,  at  least  one 
of  the  volume  -generalized  Griineisen  parameters 
must  be  negative  at  that  temperature  since  the 
model  heat  capacities  C,  and  C,  are  always  posi- 
tive. Barron2’4  and  Blackman"  have  suggested 
theoretical  models  which  for  some  volume  gener- 
alized Ouneisen  parameters,  associated  with  cer- 
tain transverse-vibrational  modes,  do  become 
negative.  These  transverse  modes  necessarily  in- 
volve elastic  constants  other  than  C„  and  Cul. 

C.  Cryogenic  measurement!  of  C|( , 

The  behavior  of  C„,  as  a function  of  temperature 
can  be  determined  directly  from  that  of  0.  Figure 
4 is  the  resulting  plot  of  Clu  versus  temperature 
for  the  four  samples.  These  curves  show  greater 
scatter  in  the  data  points  than  was  experienced  in 
previous  experiments  with  copper25  and  german- 
ium.28’27 Further,  the  estimate  of  random  error 


due  to  temperature  changes  in  the  transducer 
bonds  is  ±3%  and  the  standard  error  is  no  more 
than  ±1%.  Hence,  the  considerable  amount  of  scat- 
ter would  seem  to  indicate  structural  changes  in 
the  samples  as  a function  of  temperature.  In  large 
part,  this  may  be  due  to  the  asymmetrical  nature 
of  the  void  surrounding  the  oxygen  atom,  which  is 
a reflection  of  the  random  nature  of  the  tetrahedral 
network.  As  the  temperature  is  lowered,  evidence 
exists24  that  the  Si-O-Si  bond  angle  changes  due  to 
the  asymmetrical  vibrations  of  the  oxygen  atoms. 
Thus,  the  tetrahedral  network  of  fused  silica  be- 
low the  quenching  temperature  is  in  a state  of 
stress.  The  resulting  states  of  strain  are  more 
accentuated  for  the  weak  impurity  bonds  and  this 
makes  them  more  efficient  electron  traps  than  the 
strong  Si-O  bonds.  This  means  that  the  effect  of 
temperature  changes  is  to  cause  variations,  espe- 
cially in  anharmonic-dependent  effects  since  it  is 
these  effects  that  are  most  sensitive  to  variation 
in  structural  symmetry  and  strain  perturbations 
introduced  by  impurity  sites. 

Although  scatter  in  the  cryogenic  measurements 
is  too  large  to  attach  a great  amount  of  signifi- 
cance to  individual  data  points,  trends  in  the  Cm- 
ver sus -temperature  data  are  apparent.  The  CIU 
values  of  Suprasil  IVl,  Suprastl  H'2,  and  Suprasil  2 
all  show  a slight  general  decrease  of  Cul  with  de- 
creasing temperature.  The  Clu  values  of  Suprasil 

1 on  the  other  hand  do  not  decrease  with  decreas- 
ing temperature.  As  with  the  room -temperature 
measurements,  Suprastl  1 (with  high  OH  content) 
seems  to  be  an  anomaly.  Again  however,  Suprasil 

2 (also  with  high  OH  content)  seems  to  show  the 
general  Cm-versus-temperature  behavior  as  Su- 
prasil W1  and  Suprasil  W2  (both  with  low  OH  con- 
tent). Following  the  explanation  given  for  the 
room-temperature  measurements,  it  is  conceiv- 
able that  the  lack  of  strong  homogeneity  in  the 
Suprasil  2 sample  gives  rise  to  large  density  fluc- 
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tuations  across  the  face  of  the  ultrasonic  beam. 

The  change  of  temperature  possibly  may  exagger- 
ate these  density  fluctuations,  thereby  giving  rise 
to  the  temperature  dependence  of  C1U  observed  for 
Suprasil  2.  The  pronounced  decrease  of  the  CM1 
value  of  Suprasil  2 (Fig.  4)  with  decreasing  tem- 
perature in  the  range  147-88  K indicates  that  a 
strong  structural  change  (and  consequently  density 
changes)  may  have  taken  place. 

D.  Conclusion 

In  conclusion,  the  present  experimental  deter- 
minations of  (3  and  C1U  as  functions  of  temperature 
allow  an  important  inference  to  be  made  regarding 
the  temperature  dependence  of  the  Griineisen  pa- 
rameter for  fused  silica.  If  the  quasiharmonic 
Debye  model  of  a solid  is  assumed  to  hold  for 
fused  silica,  then  the  relatively  weak  temperature 


dependence  of  f)  and  Cm  indicates  that  elastic  con- 
stants, other  than  Cu  and  Cln,  must  dominate  the 
strong  temperature  dependence  of  the  experimen- 
tally determined  Griineisen  parameter  for  fused 
silica  at  low  temperatures.  This  is  consistent  with 
the  suggestion  of  White  and  Birch3"  that  transverse 
vibrations  of  fused  silica  are  associated  with  the 
oxygen  atoms  of  the  tetrahedral  network.  The 
clastic  constants  associated  with  these  transverse 
modes  would  be  of  the  type  C44,  Cl44,  and  C456, 
which  do  not  enter  into  the  measurements  de- 
scribed here. 
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SOMMAIRE 


ACOUSTIQUE  NONLINEAIRE  ET  PHYSIQUE  DU  SOLIDE 

L'etude  de  l'acoustique  non  lin£aire  en  plus  de  developpements 
mathematiques  et  duplications  pratiques  pleins  de  promesses  nous 
conduit  a des  progres  essentiels  dans  l'etude  de  la  physique  du  solide. 

La  deformation  non  lineaire  d'une  onde  ultra-sonore  dans  un  solide  est 
controlee  par  un  parametre  nonlineaire  qui  est  a la  fois  fonction  des 
constantes  d'elasticite  du  troisieme  ordre  [TOE]  et  fonction  des  constantes 
d'dlasticit£ — plus  connues — du  deuxifeme  ordre  [SOE] . La  mesure  de  la 
deformation  des  ondes  ultrasonores  rend  alors  possible  l'gvaluation  des 
constantes  TOE,  les  constantes  SOE  pouvant,  elles,  etre  evaluees  d'apres 
les  vitesses  des  ondes  ul trasonores. 

Le  nombre  des  constantes  elastiques  est  determine  par  la  symetrie 
du  cristal.  La  grandeur  et  le  signe  des  constantes  TOE  sont  determines 
par  les  forces  intermoleculaires.  Si  Ton  considere  un  cristal  cubique 
dans  lequel  existent  les  forces  centrales  et  les  interactions  les  plus 
proches,  a des  temperatures  voisines  de  0 °K  on  trouve  que  C^i  = = 

2C166  et  C123  = C456  = C144  = °* 

Au  cours  des  dernieres  annges  nous  avons  mesure  les  constantes  TOE 
en  fonction  de  temperature,  descendant  jusqu'a  3 °K  dans  des  cristaux 
tels  que  le  germanium  et  le  cuivre  et  dans  de  la  silice  amorphe  fusible. 
Nous  presentons  un  resume  des  resultats  accesibles  a l'heure  actuelle 
et  nous  faisons  une  comparaison  du  comportement  des  forces  inter- 
moleculaires pour  des  cristaux  de  types  differents. 

Recemment  nous  avons  decouvert  que  dans  le  cuivre  les  forces 
centrales  et  les  interactions  les  plus  proches  semblent  prevaloir  pour 
determiner  les  constantes  TOE  pour  les  temperatures  de  40  °K  et  de  190  °K. 
Toutefois  il  n'en  est  rien  a 0 °K  contrairement  ce  a quoi  l'on  pourrait 
s'attendre. 
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ABSTRACT 

NONLINEAR  ACOUSTICS  AND  SOLID  STATE  PHYSICS 

In  addition  to  mathematical  advances  and  very  promising  practical 
applications,  the  study  of  Nonlinear  Acoustics  now  is  leading  to  funda- 
mental advances  in  Solid  State  Physics.  Nonlinear  distortion  of  an 
ultrasonic  wave  in  a solid  is  controlled  by  a nonlinearity  parameter 
which  is  a function  of  the  third-order  elastic  (TOE)  constants  as  well 
as  the  second-order  elastic  (SOE)  (the  usual)  elastic  constants. 

Measurement  of  the  waveform  distortion,  then,  makes  possible  the 
evaluation  of  the  TOE  constants,  because  the  SOE  constants  can  be 
evaluated  from  ultrasonic  wave  velocities. 

The  number  of  elastic  constants  is  determined  by  the  crystal  symmetry. 
The  magnitude  and  sign  of  the  TOE  constants  are  determined  by  inter- 
molecular  forces.  In  a cubic  crystal  in  which  central  forces  and  nearest- 
neighbor  interactions  exist  one  would  find  - ZC-)-^  = ^166  anc* 

C123  = C456  = ^i 44  = 0 i n the  limit  of  0 °K. 

Over  the  past  several  years  we  have  measured  TOE  constants  as  a 
function  of  temperature  down  to  3 °K  in  such  crystals  as  germanium  and 
copper  and  in  amorphous  fused  silica.  A summary  of  presently  available 
data  is  presented,  and  a comparison  of  the  behavior  of  the  intermolecular 
forces  for  different  types  of  crystals  is  made.  Recently  we  have  found 
that  in  copper  central  forces  and  nearest-neighbor  interactions  seem  to 
predominate  in  determining  TOE  constants,  but  this  does  not  happen  near 
0 °K  as  expected.  It  seems  to  happen  near  40  °K  and  again  at  190  °K. 
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When  a sinusoidal  ultrasonic  wave  of  finite  amplitude  propagates 


I 


through  a nonlinear  medium  it  undergoes  waveform  distortion,  as  we  all 
know.  The  amount  of  this  distortion  and  its  dependence  upon  frequency 
and  propagation  distance  are  determined  by  the  nonlinearity  parameter 
of  the  propagating  medium.  For  gases  this  nonlinearity  parameter  is 
simply  y + 1 where  y is  the  ratio  of  specific  heats.  For  liquids  it  is 
B/A  + 2,  where  B/A  is  the  ratio  of  two  terms  in  the  equation  of  state. 

For  solids  the  situation  is  somewhat  more  complicated.  Not  only  are 
there  both  longitudinal  and  transverse  waves  to  consider,  but  also  the 
fact  that  solids  can  be  anisotropic  means  that  the  longitudinal  and 
transverse  waves  in  general  are  coupled,  so  that  one  cannot  define  a 
nonlinearity  parameter  for  every  possible  direction  in  a crystalline 
solid.  Our  solution  to  this  problem  has  been  to  consider  the  propagation 
of  finite  amplitude  ultrasonic  waves  only  along  the  principal 
crystallographic  directions,  and  to  restrict  our  attention  for  the 
present  to  crystals  of  cubic  symmetry.  In  the  principal  directions  pure 
mode  propagation  exists  for  the  longitudinal  wave,  even  though  the 
transverse  wave  is  always  coupled  to  a longitudinal  wave.  The  pure-mode 
longitudinal  wave  is  coupled  to  its  own  second  harmonic  in  exactly  the 
same  way  a longitudinal  wave  in  a fluid  is.  For  these  waves  a nonlinearity 
parameter  can  be  defined,  and  it  is  found  that  the  nonlinearity  parameter 
is  a function  of  the  second-order  elastic  (SOE)  constants  and  third-order 
elastic  (TOE)  constants.  The  study  of  nonlinear  distortion  of  ultrasonic 
waves  in  solids,  then,  has  produced  a means  of  measuring  the  TOE  constants 
(the  SOE  constants  can  be  determined  from  the  velocity  of  ultrasonic 
waves  in  the  sample).  Furthermore,  the  technique  can  be  adapted  to 
measurements  as  a function  of  temperature,  which  produces  results  which 
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are  of  fundamental  importance  to  solid  state  physics.  At  present  this 
is  the  only  technique  which  is  consistently  producing  such  results.  Some 
of  our  recent  results  are  quite  exciting;  however,  first  let  me  remind 
you  of  the  technique. 

In  the  principal  directions  in  a cubic  crystal  longitudinal  waves 

of  finite  amplitude  are  described  by^ 

2 2 

••  i/  3 U /oi/  , i/  \ 9U  3 U / -i  \ 

P0U  = K2  — 2 + (3K2  + K3J  9?  77  (1  > 

3d  3d 

where  and  K3  are  linear  combinations,  respectively,  of  the  SOE  and 

2 

TOE  constants,  as  shown  in  Table  I.  The  nonlinearity  parameter  is  the 
negative  of  the  ratio  of  the  nonlinear  term  to  the  linear  term  in  Eq.  1, 


namely. 


3K2  + K3 


For  an  initially  sinusoidal  disturbance  at  a = 0,  the  solution  takes 


the  form 


A?  k2  a 3 


u = A.|  sin(ka-wt)  + g cos  2(ka-wt)  + ...  (.3) 

in  which  the  second  harmonic  amplitude 

A?  k2  a3 

A2  . -J-j (4) 

contains  the  nonlinearity  parameter  3.  Our  measurements,  then,  consist 
of  absolute  measurement  of  the  second  harmonic  amplitude,  the  funda- 
mental amplitude,  the  frequency,  and  the  sample  length.  The  only  real 

problem  in  these  measurements  is  the  fact  that  the  second  harmonic 

-2  ° 

amplitude  is  so  low  (of  the  order  of  10  A at  60  MHz)  that  care  must 

be  exercised  in  getting  an  accurate  value.  A capacitive  detector  has 

worked  well  for  us,  so  that  we  have  been  able  to  determine  3 for 

3 

several  samples— not  only  at  room  temperature,  but  also  at  low  tem- 

4 

peratures— down  to  3 °K.  The  behavior  of  3 as  a function  of  temperature 


and  the  behavior  of  the  TOE  constants  as  a function  of  temperature  are 
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of  most  interest  at  present.  Let  us  contrast  their  behavior  in 
germanium,  in  which  covalent  bonding  among  the  atoms  exists,  with  that 
in  copper,  in  which  exchange  energy  resulting  from  overlapping  of  closed 
electronic  shells  seems  to  predominate  in  determining  the  elastic 
constants. 

Nonlinearity  parameters  measured  as  a function  of  temperature  for 
germanium  are  shown  in  Fig.  1;  those  for  copper  are  shown  in  Fig.  2.  In 
each  case  a smooth  curve  is  drawn  connecting  the  data  points.  As  can  be 
seen,  the  nonlinearity  parameter  is  different  in  the  three  principal 
directions  in  both  cases,  but  the  difference  in  the  magnitudes  between 
germanium  and  copper  is  not  great.  Furthermore,  the  magnitude  increases 
in  f ie  order  [100],  [111],  and  [110]  in  both  cases.  The  variation  with 
temperature  is  not  great,  the  greatest  variation  being  only  20%  for  the  J 

[110]  direction  in  germanium.  One  can  conclude,  therefore,  that  the 
nonlinearity  parameter  of  solids,  like  the  firuneisen  parameter,  should 
be  nominally  independent  of  temperature. 

There  is  much  more  information  to  be  obtained  from  these  data, 

however.  Even  though  we  cannot  isolate  all  six  of  the  TOE  constants, 

we  can  isolate  certain  combinations  of  them.  These  are  the  quantities 

2 

listed  in  Table  I as  «3.  Using  Eq.  (2),  and  known  values  of  = pV  , 
we  can  determine  the  combinations  «3  plotted  in  Figs.  3 and  4.  The  value 
of  for  the  [100]  direction  is  the  single  TOE  constant  The 

values  of  K,  in  the  other  directions  are  the  more  complicated  linear 

J * 

combinations  in  Table  I.  Upon  examination  it  is  clear  that  the  curves 
in  Figs.  3 and  4 still  are  not  the  simplest  combinations  of  TOE  constants 
available  from  the  data.  appears  in  each  expression  for  K3,  so  it 

can  be  subtracted  out.  In  this  manner  one  can  obtain  the  curves  shown 
in  Figs.  5 and  6,  which  are  the  simplest  combinations  of  TOE  constants 
available  from  our  data  on  the  nonlinearity  parameters. 
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Further  Interpretation  of  these  curves  is  facilitated  by  comparison 
with  a model.  We  will  use  the  simplest  model  available:  a model  in 
which  we  assume  that  the  atoms  interact  by  central  forces,  and  in  which 
there  are  only  nearest-neighbor  interactions.  Such  a model  has  proved 
to  be  inadequate  to  explain  the  magnitude  of  SOE  constants;  as  we  know. 

It  predicts  that  = 2C12  = 2C44,  which  is  not  satisfied  by  any  of 

the  cubic  crystals.  On  the  other  hand,  we  are  concerned  with  the  TOE 
constants  which  are  determined  by  the  nonlinear  interactions.  For  the 
TOE  constants  the  central  forces,  nearest-neighbor  interaction  model 

predicts  that  C-| -j -|  - 2C^-|2  = 2^]g6  ant^  ^123  = ^456  = ^144  = can 

use  these  expressions  to  see  how  well  those  quantities  we  can  measure 
actually  follow  the  predictions  of  the  model.  In  Figs.  5 and  6 the 
combination  C^2g  + 6C^44  + 8C^5g  should  vanish.  Fig.  5 shows  that  this 
actually  happens  at  approximately  15  °K  for  germanium.  The  fact  that 
this  combination  becomes  positive  below  this  temperature  cannot  be 
explained  at  the  moment.  Another  test  of  the  agreement  between  the 
data  and  the  central  forces,  nearest-neighbor  model  comes  from  the  fact 
that,  according  to  the  model  the  combination  C112  + 4C166  = 5/2  Clll* 

On  Figs.  5 and  6 the  value  5/2  C-j^  Is  Indicated  by  a dotted  line.  For 
germanium  the  combination  C-p2  + ^166  indeed  Is  approaching  5/2  C^-j 
near  0 #K.  This  means  that  near  0 °K  the  central  forces,  nearest-neighbor 
model  can  be  used  to  predict  quite  well  the  relationship  among  the  TOE 
constants  of  germanium.  This  Is  consistent,  because  the  model  is  expected 
to  be  most  accurate  at  0 °K,  and  less  accurate  as  the  temperature  Is 
increased. 

Let  us  now  turn  to  the  data  for  copper.  The  same  combinations  of 

/ 

TOE  constants  are  shown  In  Fig.  6.  C^-j  has  more  or  less  the  same 
behavior  for  copper  as  for  germanium,  with  only  a slightly  larger 
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variation  with  temperature.  Likewise,  there  is  reasonable  agreement 
between  + 4C.jgg  and  5/2  C^-|,  except  that  for  copper  the  agreement 
remains  essentially  the  same,  a difference  of  approximately  6%,  for  all 
temperatures.  For  copper  the  combination  C-^  + 6C^  + 8C456  1S  more 
nearly  zero  for  the  entire  temperature  range.  In  fact,  it  vanishes 
completely  at  two  temperatures:  40  °K  and  again  at  200  °K. 

Let  us  assume  for  the  moment  that  the  TOE  constants  of  copper  actually 
do  follow  the  predictions  of  the  central  forces,  nearest-neighbor  model. 

In  that  case  we  would  be  left  with  the  necessity  of  explaining  the  dips 
in  the  data  in  Fig.  6 which  occur  between  40  and  200  °K.  If  we  examine 
known  physical  mechanisms  we  find  that  dislocation  interaction,  first 

5 

reported  by  Bordoni , can  produce  in  copper  a peak  in  attenuation  which 
has  a similar  temperature  behavior.  Seeger  and  Mann^  have  pointed 
out  that  nonlinearity  is  very  large  near  a dislocation.  Therefore,  we 
are  led  to  speculate  that  we  are  viewing  for  the  first  time  the  effect 
of  Bordoni  relaxation  on  the  measured  value  of  the  TOE  constants  of 
copper.  The  samples  we  used  for  these  measurements  were  neutron 
irradiated  to  reduce  dislocation  interactions.  Thus,  if  this  is  the 
correct  interpretation  of  the  data,  we  can  say  immediately  that  our 
measurements  are  especially  sensitive  to  the  presence  of  dislocations. 
Additional  measurements  after  further  neutron  irradiation  or  after 
annealing  will  be  required  to  confirm  whether  this  is  the  correct  inter- 
pretation of  the  temperature  dependence  of  the  TOE  constants  of  copper. 

» 
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Table  1.  K2  and  for  principal  directions  of  a cubic  crystal 


Direction 

K2 

l<3 

[100] 

C11 

Clll 

[110] 

7(C11  + C]2  + 2C44) 

T^lll  + 

3C1 12  + 

12C166> 

[111] 

3<C11  * ?C12  * 4C«’ 

l(Clll  + 

6C 1 1 2 + 

1 2C1 44 

+ 24C166 

+ 2C123 

+ 16C456^ 

V 


- ■**«*..  W 
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FIGURE  CAPTIONS 


1.  Temperature  dependence  of  the  nonlinearity  parameters  of  germanium. 

2.  Temperature  dependence  of  the  nonlinearity  parameters  of  copper. 

3.  TOE  constant  combination  for  germanium  calculated  from  nonlinearity 
parameters. 

4.  TOE  constant  combination  for  copper  calculated  from  nonlinearity 
parameters. 

5.  Comparison  of  simplest  TOE  constant  combination  available  from 

germanium  data  with  predictions  of  central  forces,  nearest-neighbor 

model:  C110  + 4C.,,  = 5/2  C. ..  and  C,__  + 6C...  + 8C,c,  = 0. 
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6.  Comparison  of  simplest  copper  TOE  constant  combination  and  predictions 
of  model. 
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